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Disclaimer 

SARAH is published under the GNU library public license0, this means that you can use it for free. We 
have tested this software and its result, but we can't guarantee that this software works correctly or that 
the physical resuhs derived using this software are correct. 

We have the following requests: 

• If you find any bug, please inform us by eMaih florian.staub@physik.uni-wuerzburg.de 

• ff you have any suggestions, what is missing or can be improved, please let us know. 

• If you use this software, please cite arXiv: 0806. 0538 
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Chapter 1 

Introduction 



Supersymmetry is one of the most popular extensions of the Standard model (SM) of particle physics 
[2 [51 [31 m [SJ [S] : It solves the hierarchy problem O [S], leads to miification of the three gauge couplmgs 
[ni uni El [H] and ofFers a candidate for dark matter [13J. 

The minimal supersymmetric Standard model (MSSM) is nowadays well studied. Every event generator 
or diagram calculator can nowadays handle the MSSM out of the box. Unfortunately, there remains a lot 
of work, if you want to change the supersymmetric model, e. g. extend the particle content, add new gauge 
groups or drop some simplifying assumptions about involved parameters. In this case, it must be checked 
that the new model is free from gauge anomalies. As second step, the full Lagrangian must be derived and 
all interactions extracted. This is complicated by the fact, that the fields in gauge eigenstates have to be 
rotated to new mass eigenstates: the rotations must be incorporated, mass matrices have to be calculated 
and diagonalized. Finally, the tadpole equations are needed to find the minimum of the potential. All 
that steps are needed just for get a rough impression of the new model. If also phenomenological studies 
should be made by using one of the existing programs, model files must be created. For the embedding 
of the model in a GUT theory and getting natural values for the parameters, the Renormalization Group 
Equations (RGEs) are needed. Furthermore, often loop corrections to the masses are demanded. 
This is exactly that kind of work SARAH was written for. SARAH just needs the gauge structure, 
particle content and superpotential to produce all Information about the gauge eigenstates. As gauge 
groups, all SU (N) groups can be handled and the superfields can transform as in any arbitrary, irre- 
ducible representation of this groups. Breaking of gauge symmetries and mixings of particles can easily 
be added. Also the gauge fixing terms can be specified, and the corresponding ghost interactions are 
calculated automatically. The 2-Loop RGEs for the superpotential parameters, the gauge couplings and 
the soft breaking parameters are derived, while the self energies are calculated at one loop level. SARAH 
can write all information about the model to I^TjT;Xfiles, or creates model files for FeynArts [14] and 
CalcHep/CompHep [151 116|. which can also be used for dark matter studies using MicrOmegas |17| . 
The intention by the development of SARAH was to make it very flexible: There is a big freedom for 
the matter and gauge sector you can handle. The work with SARAH should be easy: Every information 
SARAH needs are specified in an easy to modify model file. Nevertheless, SARAH is also fast: A exist- 
ing model can be changed within minutes, and the needed time for doing all necessary calculations and 
writing a model file is normally less than 10 minutes. 

This manual is structured as follows: In the next chapter we explain the installation and the general 
setup of SARAH. In the third chapter the commands of SARAH are explained before in two foUowing 
chapters all possibilities of defining models and parameters are shown in detail. The last two chapters 
are about the possible output and the usage of numerical values. 
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Chapter 2 

Qmck Start 



2.1 Download and Installation 

SARAH is a package for Mathematica and was tested with versions 5.2 and 7.0. 
SARAH can be downloaded from 

http : / / theorie . physik . uni-wuerzburg . de/~f nstaub/ sareih . html 

The package archive contains the following directories: 

• Models: AU models you want to run with SARAH, have to be saved here. This directory contains 
aheady some models hke MSSM and NMSSM. 

• Package: AU package files. 

• LaTeX-Packages: I^TjTjKpackages, which might be needed for the output of SARAH, if there are 
not instaUed on your machine already. 

During the work, there is also the directory 

Output 

created. It wiU contain aU files written by SARAH. 

You can save this file of course everywhere on your machine, but the most convenient choice would be in 
the application directory of Mathematica. This is under Linux 

home/user/ .Mathematica/ Applications/ 

and 

Mathematica-Directory\AddOns\Applications\ 
in Windows. 

^Mathematica is a protected product by Wolfram Research 
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2.2 Run S ARAH 

After you have extracted and saved the package directory, you can load SARAH in Mathematica in 
different ways. If the package directory is in the apphcation directory of Mathematica, you can just use 

«"sarah-2 . 0/SARAH.m" 

Otherwise, you have to use the complete path 
«"/path/SARAH.m" 

After SARAH is successfully loaded, you start the calculations for a specific model with 
Start ["Modelname"] ; 

Here, Modelname is the name of a subdirectory of the Model-directory, e. g. for the minimal supersym- 
metric Standard model use 

Start ["MSSM"] ; 

It is also possible to define several submodels of one model, e. g. the MSSM in CKM basis, the MSSM 
with flavor violation or without flavor violation. These models must be saved in a subdirectory of the 
corresponding model directory and can be chosen via 

Start ["Modelname" , "Submodelname"] ; 

For example, the MSSM without flavor violation is loaded via 
Start ["MSSM", "NoFV"] ; 

2.3 What happens automatically 

After you have used Start [modelf ile] , SARAH first checks, if the selected model is free of gauge anoma- 
lies (see appendix [D]) . If there is a gauge anomaly, a warning is printed. After that, SARAH begins with 
its main work and calculates the different parts of the Lagrangian for the gauge eigenstates from the 
given gauge groups and superpotential. 

The next steps are to accomplish all necessary rotations and redefinition of fields: If a gauge symmetry 
is broken, the fields responsible for the symmetry breaking are getting a vacuum expectation value and 
the gauge fields are rotated. Afterwards the matter particles are rotated to the new mass eigenbasis and 
the tadpole equations are derived. This steps can be repeated, if more rotations or symmetry breakings 
are necessary. 

During this evaluation some more things might done automatically: Particles are integrated out to get 
an effective theory, the ghost interactions are derived from the gauge fixing part of the Lagrangian or 
the parametrization of defined interactions is changed. At the end, SARAH calculates for all appearing 
particles in the model the masses at tree level and splits the Lagrangian in different pieces to increase 
the speed of following calculations. Additionally, if numerical values for all parameters are provided, e. g. 
in the input files or in a Les Houches spectrum file [H], SARAH calculates the eigenvalues of the mass 
matrices and the rotation matrices. 

When the evaluation of the model is finished, the user has the choice if he wants to work directly with 
SARAH, e. g. playing with mass matrices or calculating vertices, or if complete model files should be 
generated by SARAH. 
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2.4 Commands 

The most important commands to work with SARAH are: 

• Vertex [Fields , Options] : Calculates a vertex for given fields 

• MassMatrices [States] : Shows all mass matrices for givcn cigenstates States 

• TadpoleEquations [States] : Shows ah tadpole equatioiis for given eigenstates States 

• CalcRGEs [Options] : Calculates the RGEs 

• CalcLoopCorrections [Options] : Calculates one loop corrections for givcn cigenstates States 

• ModelOutput [States, Options]: Create output defined by options for given eigenstates States. 

• MakeVertexList [States, Options]: Calculates all vertices for given eigenstates States 

• MakeTeX [Options] ; Writes M^files 

• MakeCHeop [Options] : Write CalcHep/CompHep model files 

• MeikeFeynArts: Write FeynArts model file 
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Particie Content, Gauge Sector and Superpotential 
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Check Gauge Anomaiies 



Basic Lagrangian 



Integrate out Gauge Eigenstates 



Add Ghost Interactions 



Start Set of Rotations 
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Redefine Couplings 
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Lagrangian of Gauge Eigenstates 




Mixings Gauge Sector 
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Mixings Mass Sector 
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Figure 2.1: Setup of SARAH 
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Chapter 3 

Working with SARAH 

3.1 Definition of Particles 

Before we explain the different functions of SARAH, we have to clarify the nomenclature. Ali gauge 
eigenstates are named as follows: 

ParticleType <> Basis <> [Indices] 

Type Here, ParticleType is just one letter and indicates the type of a field. The conveiition is as 
follows: 

• F for fermionic compoiient of chiral superfield 

• S for scalar component of chiral superfield 

• f for fermionic component of vector superfield 

• V for bosonic component of vector superfield 

• g for ghost field 

In addition, there are two types of auxiliary fields. This auxiliary field are not related in any way to 
the auxiliary components of the superfields in SUSY theories, but they are needed only for writing a 
CalcHep/CompHep file (see l6.3.3)) . The type indicating letters are: 

• A for an auxiliary scalar 

• a for an auxiliary vector boson 

After rotating fields in new mass eigenstates, there are no any longer constraints 
on the names of fermions and scalars, but vector bosons must still begin with 
V, because the corresponding ghosts are automatically added! 



Basis The Basis of a particle is the name of the underlying superfield. In the MSSM for example this 
can be dL for the superfield of the left handed down quarks and squarks or G for the vector superfield 
transforming under the strong interaction. 
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Indices Indices is the list of the indices the particle carries. There are three difFerent kinds of indices: 

• generation: For all particles, which appear in more than one generation 

• lorentz : For all particles carrying a Lorentz iiidex 

• charge : For all particles charged under a non-abelian gauge group, if tliis indices are not expanded 
(I5.1.1.U.1I) . 

Examples To clarify the above definitions, here some exainples: 

a) VB [{lorentz}] : B-Boson with one Lorentz index 

b) f B; Bino 

c) SHdO; Neutral down Higgs 

d) FHdO; Neutral down Higgsino 

e) VG [{generation, Lorentz}] : Gluon with one generation and one Lorentz index 

f) f G [{generation}] : Gluino with one generation index 

g) gG [{generation}] : Gluon ghost with one generation index 

h) SdL [{generation, color}] : Left handed d-squark with one generation and one color index 
i) FdL [{generation, color}] : Left handed d-quark with one generation and one color index 
j) hh [{generation}] : Neutral, CP-even Higgs (light and heavy Higgs) after EWSB 

For a list of all particles defined in the models delivered with SARAH see appendix [e1 

3.1.1 Antiparticles 

There are two function to assign antiparticles: Depending on the type of the particle conj or bar are 
used. 

• Scalar, vector boson and Weyl spinor: conj , e. g. conj [SdR] or conj [VWm] . 

• Dirac fermion and ghost: bar, e. g bar [Fd] or bar [gG] . 

SARAH checks every time, if a particle is real or complex or a majorana fermion and simplihes the 
expressions by using 

conj [RP] : = RP 
bar [MF] := MF 

for a real particle RP or real parameter, and for a majorana fermion MF. The name of all majorana fermions 
of the current model is saved in the list MajoranaPart, and all real parameters and particles are listed 
in realVar. 

Some words about the necessity of conj : The function conj is very similar to the existing function 
Conjugate of Mathematica. But at least with Mathematica 5.2 there are some problems concerning this 
function: It is not possible to calculate the derivative with respect to a complex conjugate variable, and 
Conjugate is a pure numerical functions and therefore in some cases too slow in handling big analytical 
expressions. This was improved in Mathematica 7.0, but we want to stay compatible also with 5.2. 



14 



3.2. DIFFERENT EIGENSTATES OF ONE MODEL 



3.2 Different Eigenstates of one Model 

While calculating the Lagrangiaii of a model, SARAH saves the information of the different eigenstates 
of the model. The name of the eigenstates can be defined by the user in the model file, see (|5.3|) . For the 
Standard model file this are 

• GaugeES: Gauge eigenstates without any rotation 

• EWSB: Eigenstates after electroweak symmetry breaking. 

• SCKM: Eigenstates in Super-CKM basis 

• TEMP: Auxiliary eigenstates with no physical meaning 

3.3 Model Information 

There is a lot of information automatically calculated and saved by SARAH. In this section is shown, 
how this information can be accessed and used. 

3.3.1 Particle Content 

To get an overview of all particles of the different eigenstates, use 
Particles ["Name of Eigenstates"] 

e.g. Particles [GaugeES] or Particles [EWSB] . 

This returns a list with the following information 

• Name of the particle 

• Type of the particle (F for fermion, S for scalar, G for ghosts, A for auxiliary field) 

• Number of first generation (can be different from 1 in effective theories) 

• Number of last generation 

• Indices of the particle 

For fermions only the Weyl spinors and not the Dirac spinors are shown. 

Example A typical entry in Particles [GaugeES] looks like 
{FdL, 1, 3, F, {generation, color}} 

3.3.2 Masses of Particles 

SARAH automatically calculates the tree level masses of all particles. This information is saved in 
Masses ["Name of Eigenstates"] 

It returns are a list with replacements rules. First a statement Mass [] with the name of the particle is 
given followed by the value. There are three possibilities: 
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• Expression: If a expression is used for a mass, SARAH has calculated the tree level mass depending 
on other parameters of the model. 

• MassGiven: This means that a numerical value for the mass was given in the particle file fsee l4.3p . 

• MassRead: This means that the value of the mass was read in from a Les Houches file. 
More Information of defining masses is given in section 14.31 

Example The mass of the Z-Boson (VZ) after EWSB is saved in Masses [EWSB] . The corresponding 
entry is 

Mass[VZ] -> ((vd-2 + vu-2) * (g2*Cos [ThetaW] + gl*Sin [ThetaW] ) "2/4 

3.3.3 Mass Matrices 

SARAH calculates automatically the mass matrices before rotating the fields in the new eigenstates and 
saves the Information in arrays. The basis of the rotations can be seen by using 

MixBasis ["Name of Eigenstates"] 

The matrices itself are saved in two arrays: 
MassMatrices ["Name of Eigenstates"] 
And 

MassMatricesFull ["Name of Eigenstates"] 

The difference between this two arrays is that in the first one, the difFerent generations are written as 
indices, while in the second on the generation indices are explicitly inserted. This means, in the first case 
the basis for the mass matrix in the down squark sector is just 

( SdL [{gn , cn}] , SdL [{gm , cm}] ) 

while in the second case the basis vector is 

(SdL[{l,cnl}] ,SdL[{2,cn2}] , SdL [{3 , cn3}] , SdL [{1 , cml>] ,SdL[{2,cm2}] , SdL [{3 , cm2}] ) 

Example Let's have a look at the down-squark sector after EWSB. 

MixBasis [EWSB] [[!]] 

returns the basis of the mass matrix: 

{{SdL, SdR}, {conj [SdL] , conj[SdR]}} 

The (1,1) element of the mass matrix is saved in 

MassMatricesFull [EWSB] [[1,1,1]] 

and looks like 

-(gl-2*vd-2)/24 - (g2-2+vd-2)/8 + (gl-2*vu-2) /24 + 
(g2-2*vu-2)/8 + mq2[l, 1] + (vd-2*Yd[l, l]-2)/2 
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3.3.4 Tadpole Equations 

Gauge symmetries are broken, if particles receive an Vacuum Expectation Value (VEV). This VEVs 
should minimize the potential and therefore fulfill the equations 

dV 

j- = (3.1) 

V is the Potential and Vi are the difFerent VEVs. This are the so called tadpole equations, which are 
cubic equations in the VEVs. SARAH saves the left hand side of (|3.ip in arrays named 

TadpoleEquations ["Name of Eigenstates"] 

The order ol the tadpole equations in this array corresponds to the order of the definition of VEVs in the 
model file, see (|5.3.1.ip . 

Example The tadpole equation for Vd after EWSB is saved in 
TadpoleEquations [EWSB] [[1]] 

and reads 

mHd2*vd + (gl~2*vd~3) /8 + (g2~2*vd~3) /8 - (gl~2*vd*vu~2) /8 - 
(g2~2*vd*vu~2)/8 + vd*\[Mu]~2 - vu*A[\[Mu]] 

3.3.5 Parts of the Lagrangian 

The most important, but also most coniplex piece of a model is the full Lagrangian. SARAH calculates 
the full Lagrangian from the superpotential and the gauge sector by using the method explained in 
appendix|Al The final results, the Lagrangians for the different eigenstates are saved as 

Lagrangian ["Name of Eigenstates"] 

For a realistic SUSY model, the Lagrangian is generally very lengthy. Therefore SARAH splits it in 
difFerent parts in order to speed up some calculations. This splitting might be also helpful for analyzing 
the structure of interactions involving different kinds of fields. An overview of all names for the difFerent 
parts oF the Lagrangian is given in IA. 61 

Also the results of the different steps during then calculation of the Lagrangian in gauge eigenstates (e. g. 
F-Terms, D-Terms, kinetic parts) are saved, so it is possible to have a detailed look at specific parts of 
the Lagrangian. The names of the parts are also given in appendix 1X1 

3.3.6 Write all information about Particles and Parameters in external File 

Ali information about particles and parameters of the considered model can by written in two files 
Particlelnf o .m and Parameterinf o, which are saved in 

. ./SARAH/Output/"Model Name"/ 

by just using 

ExportModelInf ormation; 

The file Particlelnf o .m contains the following information for all eigenstates of the model 
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• R-Parity 

• PDG 

• Wr^name 

• Output name 

• FeynArts number 

• Type 

• Self-conjugated or not 

• Generatioiis 

• Indices 

while Pareuneterlnf o .m contains the Information about all parameters 

• Dependence on other parameters 

• Real or complex 

• Numerical value 

• Les Houches input 

• I^TEKname 

3.4 Calculating Vertices 

One main function of SARAH is to calculate the vertices for a model. In contrast to the most other 
calculations, vertices normally are not calculated automatically: It can last several minutcs to calculate 
all vertices of a model and if one only wants to check some interactions. this calculations are not necessary. 
Of course, it is also possible to tell SARAH that all vertices should automatically calculated, see section 
16.41 In this section we want to focus on calculating vertices 'by hand'. 

Calculating vertices is done with 
Vertex [ParticleList ,Qptions] 

ParticleList is an list containing the involved fields. This list can consist of up to 6 particles, if you 

analyzing a theory with effective operators (see I5.4p . 

The following Options are supported by the Vertex command: 

• Eigenstates, Value: Name of Eigenstates, Default: EWSB 
Fixes the considered eigenstates 

• UseDependences: Value True or False, Default: False 

Optional relations between the parameters (see section [4. 2 p will be used, if UseDependences is set 
to True. 

The Output of Vertex is an array: 
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■[{PeirticleList},{{Coefficient 1, Lorentz l},{Coef f icient 2, Lorentz 2},...} 

First, the list of the involved particles is given and thc indiccs aro inscrtcd. The second part consists of 
the value of the vertex and can be also a hst, if different Lorentz structures are possible. In the part 
independent of Lorentz indices can appear 

• Delta [a, b]: Kronecker Delta (5«^ 

• Lam [t , a , b] : Gell-Mann Matrix A^^ 

• LambdaProd [x , y] [a, b] : Matrix product of two Gell-Mann Matrices (A^A^)^^ 

• Sig[t,a,b]: Pauli Matrix 

• SigmaProd[x,y] [a, b] : Matrix product of two Pauli Matrices (o'^cr^)^^ 

• fSU3[i,j,k]: Structure constants of 5C/(3): f^'' 

• fSU2[i,j,k]: Structure constants of S'C/(2): e'^'' 

• FST[SU[N]] [i,j ,k]: Structure constants of 5C/(7V): f j}'' 

• Couplings: gl, g2, g3, Ye [a,b] , Yd[a,b] , Yu[a,b] , . . . 

• Mixing matrices: ZD[a,b], ... 
While the Lorentz part can consist of 

• gamma [lor] : Gamma matrix 

• g[lorl,lor2] : Metrik g^,^ 

• Mom [particle , lor] : Momentum pp of particle P 

• PL, PR: Polarization operators Pl = Pr = ^-^^ 

• 1: If the vertex is a Lorentz scalar. 

• LorentzProduct [_,_]: A non commutative product of terms transforming under lorentz group 

Examples Some examples to clarify the usage and output of Vertex: 
a) Scalar - Scalar - Vector 

Vertex[{hh,Ah,Z}] 

leads to the vertex of Higgs - Pseudo scalar Higgs - Z: 

{{hh[{gtl}], Ah[{gt2}], VZ[{lt3>]>, 
{((MA[gt2, l]*MH[gtl, 1] - MA[gt2, 2]*MH[gtl, 2])* 

(g2*Cos [ThetaW] + gl*Sin[ThetaW] ))/2, 
Mom[Ah[-[gt2}] , ItS] - Moin[lili[{gtl}] ,lt3] }} 

We see that the output is divided in two parts: First, the involved particles are given, second, the 
value of the vertex is given. This second part is again split in two parts: The eoefBcient independent 
of Lorentz indices and the part defining the transformation under the Lorentz group. 
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b) Fermion - Fermion - Scalar 
Vertex [{bar [Fd] .Fd.hli}] 

is the interaction between d-quarks and Higgs: 

{{bar[Fd[{gtl, ctl}]], Fd[{gt2, ct2}] , hh[{gt3}]}, 
{((-I)*Delta[ctl, ct2]*Delta[gtl, gt2]*MH[gt3, 2]*Yd[gt2, gtl] ) /Sqrt [2] , PL>, 
{((-I)*Delta[ctl, ct2]*Delta[gtl, gt2]*MH[gt3, 2]*Yd[gtl, gt2] ) /Sqrt [2] , PR» 

Obviously, there are three parts: One for the involved particles and two for the different Lorentz 
structure. 

c) Weyl Fermions - Vector Boson 

It is also possible to calculate the vertices for Weyl and not for Dirac fermions: 

Vertex[{fB, FdL, conj[SdL]}, Eigenstates -> GaugeES] 

gives 

{{fB, FdL[{gt2, ct2}] , conj [SdL[{gt3, ct3}]]}, 
{((-I/3)*gl*Delta[ct2, ct3] *Delta[gt2, gt3] ) /Sqrt [2] , 1}} 

d) Ghost interactions 

Vertex [{bar [gWm] .gWm.hh}] 
leads to 

{{bar[gWm], gWm, hh[{gt3}]}, 
{ (-1/2) *g2*Mass[VWin]*(Cos[\ [Beta] ]*MH[gt3, 1] - MH[gt3, 2] *Sin[\ [Beta] ] ) , 1}} 

e) Using of Dependences: 
With 

Vertex [{conj [Se] , Se , VP} , UseDependences -> True] 

gi and 52 are replaced by the electric charge e: 
{{conj [Se [{gtl}]], Se[{gt2}], VP[{lt3}]}, 

{(-I)*e*Delta[gtl, gt2] , -Mom [conj [Se [{gtl}]] , lt3] + Moin[Se [{gt2}] , lt3]}} 

f) Fixing the generations 

It is possible to fix the indices of the particles already by using the Vertex command: 

Vertex[{hh[{l}] , hh[{l}] , Ah[{2}] , Ah[{2}]}] 
leads to 

{{hh[{l}], hh[{l}], Ah [{2}], Ah [{2}]}, 
{(-I/4)*(gl-2 + g2-2) *Cos[2*\ [Alpha] ]*Cos[2*\ [Beta] ] , 1}} 
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Obviously, the given definition of the mixing inatrix is already inserted. If you have instead fixed 
the indices by a replacement, the definition of the niixing matrix wouldn't be used: 

Vertex[{hh, hh, Ah, Ah}] /. {gtl->l, gt2->l,gt3->2, gt3->2} 

returns 

{{hh[{l}], hh[{l}], Ah [{2}], Ah[{gt4}]}, 

{(-I/4)*(gl-2 + g2-2)*(conj [ZA[2, 1] ] *conj [ZA [gt4, 1]] - 

conj[ZA[2, 2]]*conj [ZA[gt4, 2] ] ) * (conj [ZH [1 , l]]-2 - conj[ZH[l, 2]] "2), 1» 

But 

Vertex[{hh, hh, Ah, Ah}] /. {gtl->l, gt2->l,gt3->2, gt3->2} / . subAlways 

leads to the former expression. 
g) Effective Operators 

In effective theories also interactions between two fermions and two scalars are possible. As example 
an effective vertex for a model in which the gluino was integrated out: 

Vertex[{Fd, Fd, conj [Sd] , conj[Sd]}] 

Returns 

{{Fd[{gtl, ctl}] , Fd[{gt2, ct2}] , conj [Sd[-[gt3, ct3}]], conj [Sd[{gt4, ct4}]]}, 
{-(g3~2*(sum[jl, 1, 8, (Lam[jl, ct3, ct2] *Lam [j 1 , ct4, ctl] ) /Mass [f G] [j 1] ] * 
ZD[gt3, gt2]*ZD[gt4, gtl] + 
sum[jl, 1, 8, (Lam[jl, ct3, ctl] *Lam [j 1 , ct4, ct2] ) /Mass [f G] [j 1] ] * 
ZD[gt3, gtl]*ZD[gt4, gt2])), 
LorentzProduct [PL, PL]}, {O, LorentzProduct [PR, PL]}, 
{g3~2*(suni[jl, 1, 8, (Lam[jl, ct2, ct3]*Lam[jl, ct4, ctl] )/Mass [f G] [j 1] ] * 
ZD[gt3, 3 + gt2]*ZD[gt4, gtl] + 

sum[jl, 1, 8, (Lam[jl, ct2, ct4]*Lam[jl, ct3, ctl] )/Mass [f G] [j 1] ] * 

ZD[gt3, gtl]*ZD[gt4, 3 + gt2]), 
LorentzProduct [PL, PR]}, {O, LorentzProduct [PR, PR]}, 
{O, LorentzProduct [gcUimia,PL, PL] } , {O, LorentzProduct [gamma, PR, PL]}, 
{O, LorentzProduct [gcUfflna, PL, PR]}, {O, LorentzProduct [gamma, PR, PR]}} 

Obviously, SARAH checks the eight possible operators (4 different combination of polarization 
operators with and without a 7 matrix) and returns the result for each operator. 



3.5 Renormalization Group Equations 

SARAH calculates the renormalization group equations (RGEs) for the parameters of the superpotential, 
the soft breaking terms and the gauge couplings at one and two loop level. This is done by using the 
generic formulas of |31j . 

The calculation is started via 
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CalcRGEs [Options] 

The different options are 

• TwoLoop, values: True or False, Default: True 
If also the two loop RGEs should be calculated. 

• ReadLists, values: True or False, Default: False 

If the RGEs have already be calculated, the results are saved in the output directory. This results 
can be read from the files instead of doing the complete calculation again. 

The results will be stored in the foUowing arrays: 

• Gij: Anomalous dimensions 

• BetaYijk: Trilinear superpotential parameters 

• BetaMuij: Bilinear superpotential parameters 

• BetaAijk: Trilinear soft breaking parameters 

• BetaBi j : Bilinear soft breaking parameters 

• Betain2ij: Scalar soft breaking masses 

• BetaMi: Gaugino soft breaking masses 

• BetaGauge: Gauge couplings 

These arrays are also saved in the directory 
. . /SARAH/Output/"ModelName"/RGE 

Ali entries of this arrays are three dimensional: The first entry is the name of the parameter, the second 

the one loop /3-function and the third one the two loop /3-function. 

Generally, the results contain sums over the generation indices of the particles in the loop. But SARAH 
always tries to write them as matrix multiplications, in order to simplify the expressions. Therefore, new 
symbols are introduced: 

• MatMul [A , B , C , . . .] [i, j]: {ABC . . Matrix multiplication, also used for vector-matrix and 
vector-vector multiplication. 

• trace [A,B,C, . . .] : Tr{ABC .. .). Trace of a matrix or product of matrices. 

• Adj [M] : . Adjoint of a matrix 

• T[M] : M'^ . Transposed of a matrix 

To difFcr bctwccn generation and other indices during the calculation, Kronecker [i, j] is used for gen- 
eration indices instead of Delta [i, j] , but has the same meaning. 
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Example The Yukawa couplings of the MSSM are saved in BetaYijk. The first entries consists of 

BetaYi jk [[1,1]]: Ye [il , 12] 

The corresponding one loop /3-function is 

BetaYijk [[1,2]] : 

(-9*gl-2*Ye[il, 12] ) /5 - 3*g2-2*Ye [il , i2] + 3*trace[Yd, Adj [Yd] ] *Ye [il , 12] + 
trace[Ye, Adj [Ye] ] *Ye [11 , 12] + 3*MatMul[Ye, Adj [Ye] , Ye][ll, 12] 

The two loop /3-function of the electron Yukawa couphng is saved in BetaYijk [ [1 , 3] ] , but we skip it 
here because of its length. 

The results for the scalar masses are simphfied by using abbreviations for often appearing traces. The 
definition of the traces are saved in the array TraceAbbr. 

Example The soft breaking mass of the electron is the first entry of Betajn21j, and the one loop 
/3-function is 

(_24*gl~2*MassB*conj [MassB]*Kronecker[ll, 12]) /5 + 2*gl~2*Kronecker [11 , 12]+Trl[l] + 
4*inHd2*MatMul [Ye , Adj[Ye]][ll, 12] + 4*MatMul [A [Ye] , Adj [A [Ye] ] ] [11 , 12] + 
2*MatMul [me2 , Ye, Adj[Ye]][ll, 12] + 4*MatMul[Ye, ml2, Adj[Ye]][ll, 12] + 
2*MatMul[Ye, Adj [Ye] , me2][ll, 12] 

The definition of the element Tri [1] is saved in TraceAbbr [[1,1]]: 

{Trl[l], -mHd2 + mHu2 + trace [md2] + trace [ine2] - trace [ml2] + trace [mq2] - 2*trace [mu2] } 

3.6 Loop Corrections 

SARAH calculates the analytical expressions for the one loop corrections to the tadpoles and the self 
energy of all particles. This calculations are performed in the t'Hooft gauge and are started via 

CalcLoopCorrectlons ["Ncune Of Elgenstates"] 

If the vertices for the given set of eigenstates were not calculated before, this is done before the calculation 
of the loop contributions begins. 

Using the conventions of [19J, the results will contain the foUowing Passarino Veltman integrals 

• A0[m-2]: Ao{m'^) 

• BO [p~2,ml~2,m2~2] : B o {p'^, ml, ml) 

• Bl[p-2,ml-2,m2-2]: Bi{p'^,ml,ml) 

• B22[p-2,ml-2, 012-2] : B22{p^,mlml) = B22{p^,ml,ml) - lAoimj) - ^Aoiml) 

• F0[p-2,inl-2,m2-2]: Ao{ml) - 2Ao{ml) - (2p2 +2m\ - ml)Bo{p'^ ,ml,ml) 

• GO [p-2,ml-2,m2-2] : {p^ - m\ ~ ml)Bf){p'^ ,m\,ml) - AQ{m\) - Af){m\) 

• H0[p-2,ml-2,m2-2]: 4B22(^|^ "i^, ^1) + G'o(p^ m?, m|) 
while the needed couplings are named: 
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• Cp[pl,p2,p3] and Cp[pl,p2,p3,p4] for non-chiral, three and four point interactions involving the 
particles pl - p4. 

• Cp[pl,p2,p3] [PL] and Cp[pl,p2,p3] [PR] for chiral, three-point interactions involving the fields 
pl - p3. 

The result for the loop correction is generated twice: 

• As sum of the different contribution hke it would appear in a textbook 

• As Hst of the different contribution. For each contribution an entry Uke this exists: 

{Particles, Vertices, Type.Charge Factor, Sjnmnetry Factor} 

— Particles: The particles in the loop. Always two a given, also for the correction of tadpoles 
or correction of propagators due to four point interactions. For this cases, both fields are of 

course the samo ono. 

— Vertices: The needed Vertex for the correction is given. When calculating the corrections to 

mass matrices, particles with 'U' at the first position do appear (e. g. Cp [conj [USd [{gOl}] ] , Sd [{gll}] , VZ] ). 

This mcans, that it is an external ficld, which must not bc rotated ('U' for unrotated). There- 
fore, in the vertex has the niixing matrix corresponding to this vertex replaced by the identity 
matrix. 

— Cheirge Factor: If several gauge charges of one particle are allowed in the loop, this factor will 
be unequal to one. For the MSSM only the a factor of 3 can appear because of the different 

colors. 

— Symmetry Factor: If the particles in the loop indistinguishable, the weight of the contribution 
is only half of the casc of distinguishable particles. The absolute value of the factor depends 
on the type of the diagram. 

The Information about the loop correction are also saved in the directory 
. ./SARAH/Output/"ModelNaine"/"Naine of Eigenstates'VLoop 

3.6.1 One Loop Tadpoles 

The complete results as sums of the different contributions are saved the two dimensional array 
OneLoopTadpoles ["Name of Eigenstates"] 

The first column givcs the name of the corresponding VEV, the second entry the one loop correction. 
A list of the different contributions, including symmetry and charge factors, is 

LoopCorrectionTadpoles ["Name of Eigenstates"]; 

Example The first correction of the tadpoles are saved in 
LoopCorrectionTadpolees [EWSB] [[1]] ; 

and read 

{bar[Cha], Cha, Cp [Uhh [{gOl}] , bar[Cha[{gIl>]] , Cha[{gI2}]], FFS, 1, 1/2} 
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3.6.2 One Loop Self Energies 

The results are saved in the following two dimensional array 

OneLoopSelf Energy ["Name of Eigenstates"] 

The first column gives the name of the particle, the entry in the second column depends on the type of 
the field 

• Scalars: One loop self energy n(p^) 

• Fermions: One loop self energies for the difFerent polarizations (S]^(p'^),I]^(p^), S'^(p^)) 

• Vector bosons: One loop, transversal self energy 11^ (p^) 

Also lists of the different contributions do exist. One loop correetion for rotated states are saved in 
LoopCorrectionMassMatrices ["Name Of States"] 

and 

LoopCorrectionUninixed["NameuOf uStates"] 
for unrotated states. 

Example 

a) Corrections to mass matrices 

The correetion to the down squark mass matrix due to pseudo scalar Higgs is saved in 

LoopCorrectionMassMatrices [EWSB] [[1, 11]] 

and reads 

■CAh, Ah, Cp[conj[USd[{g01}]], USd[{g02}] , Ah[{gll}] , Ah[{gI2}]], SSSS, 1, 1/2} 

First, the fields in the loop are given. Afterwards, the needed vertex with indices (gOi for outer 
particlcs. gli for internal particles). The U assigns, that the mixing matrix of the squarks has to 
be replaced by the identity matrix. The type of the correetion is SSSS, i. e. a scalar-four point 
interaction, the color factor is 1 and the symmetry factor is 5, because the pseudo scalars Higgs 
are real fields. 

b) Correetion to masses The corrections to the Z mass are saved in 
LoopCorrect ionUninixed [EWSB] [ [10] ] 

An arbitrary entry looks like 

{bar[Fd], Fd, Cp[VZ, bar [Fd [{gll}] ] , Fd[{gI2>]], FFV, 3, 1/2} 



25 



Working with SARAH 



3.7 WF^- and Modelfile-Output 

With 

ModelOutput [Eigenstates , Options] 

you can produce I^TEX-files and model files for FeynArts and CompHep/CalcHep. Here Eigenstates 
specifies the eigenstates of your model, for which the output should be generated, e. g. GaugeES or EWSB. 
The Options are the following 

• WriteTeX, values: True or False, Default: False 

If a I^T[5Xfile containing the particles, mass matrices, tadpole equations and all vertices of the model 
should be written. 

• WriteFeynArts, values: True or False, Default: False 
If a model file for FeynArts should be written. 

• WriteCHep, values: True or False, Default: False 

If a model file for CompHep/CalcHep should be written. 

• ef f ectiveOperators, Values: True or False, Default: False 

If for efFective theories also higher dimensional operators should be calculated. By default this are 
only four point interactions. 

• SixParticleInteractions, Values: True or False, Default: False 

If for effective theores also the six point interactions should be calculated. 

• FeynmanDiagrcuns, Values: True or False, Default: True 

If in the I^T^r^Cfile Feynman diagrams for each vertex should be drawn. 

• ReadLists, Values: True or False, Default: False 

If before the command ModelOutput was used, lists with all vertices do exist and were saved in the 
output directory. So it is possible to read this lists and save the time for the calculation. 

• IncludeRGEs, Values: True or False, Default: True 
If the RGEs should be calculated. 

• TwoLoopRGEs, Values: True or False, Default: True 

If the two loop RGEs should be calculated. (IncludeRGEs must be True) 

• IncludeLoopCorrections, Values: True or False, Default: True 

If the one loop corrections to the self energy and the tadpoles should be calculated. 

The generated output will be saved in 

. ./SARAH/Output/"ModelName"/"Name of Eigenstates"/TeX/" 

. ./SARAH/Output/"ModelName"/"Name of Eigenstates "/FeynArts/" 

. ./SARAH/Output/"ModelName"/"Name of Eigenstates "/CHep/" 

More details about the output is given in chapter |51 
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AU information of the different models are saved in three different files, which have to be in one directory 
. . ./SARAH/Models/"ModelName"/ 



The directory name is equal to the name of the model and must contain a model file with the same 
name! 



The three files are: one model file with the same name as the directory (ModelNeune.m), a file containing 
additional information about the particles of the model (particles .m) and a file containing additional 
information about the parameters of the model (parameters .m). Only tlic first file is really necessary for 
calculating the Lagrangian and gct a first impression of a model, but for defining properties of parameters 
and particles and for producing an appropriate output the other two files are needed. 



4.1 The Model File 



The model file is the heart of SARAH: The complete model is specified by the model file. Therefore we 
will explain the general structure of the model file here, and have a look at the different functions and 
its physical meaning in detail in the next chapter. 



4.1.1 Description 

The model file contains different parts: First the gauge structure and the particle content are given, and 
the matter interactions are defined by the supcrpotential. This are general information needed for all 

eigenstates of tho model and must always be apparent. 

This part is followed by the definition of the names for all eigenstates. For these eigenstates several 
properties can be defined: decomposition of scalars in scalar, pseudoscalar and VEV, rotations in the 
matter and gauge scctor and the corresponding gauge fixing terms. New couplings can be added and 
existing couplings can be changed by hand. 

Afterwards, the particles are states, which should be integrated out or deleted. At the end, the Dirac 
spinors have to be built out of Weyl spinors, a spectrum file can be defined and a choice for automatically 
output can be made. 
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4.1.2 Schematic Structure 

The model file is structured as follows: 

General Information 

a) Gauge Structure of the model given by the vector superfields, see lS.l.ll 
Gauge[[l]] = {. . . 

b) Matter Conteut given by the chiral superfields see 15.1.21 
Fields[[l]] = { . . . 

c) Superpotential: Tri- and bilinear Terms, see 15.21 
SuperPotential ■[ ... 

Names of Eigenstates 

Names for the different eigenstates, see 15.31 
NameOf States = { . . . 

Defining Properties of the different Eigenstates 

a) Vacuum Expectation Values, see l5.3.lTT] 
DEFINITION["Name of Eigenstates"] [VEVs] = { ... 

b) Rotation of Gauge Bosons, see 15.3.21 
DEFINITION["Name of Eigenstates"] [GaugeSector] = i ... 

c) Rotation of fields, see 15.3.31 

DEFINITION["Name of Eigenstates"] [MatterSector] = { ... 

d) Gauge Fixing part, see IA. 51 

DEFINITION["Name of Eigenstates"] [GaugeFixing] = i ... 

e) Additional couplings or redefinition of existing couplings, see l5.3.5T2l 
DEFINITION["Name of Eigenstates"] [Additional] = { ... 

Additional, general Information 



a) Define Dirac Spinors, see 15.5 
dirac[[l]] = { . . . 
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b) Iiitegrating out or Deleting of Particles, see 15.41 

IntegrateOut[[l]] = { ... 
DeleteParticles[[l]] = { ... 

c) Automatized Output, see 16.41 
makeOutput = { . . . 

d) Define spectrum file, see l7.ll 
SpectrumFile = ... 



4.2 Parameter File 
4.2.1 General 

The information of the parameter file are needed for some in- and output routines. Also, you can apply 
assumptions to the parameters. The parameter file consists of a list called ParsuneterDef inition. This 
is an array with two columns: First column gives the name of the parameter, second column defines the 
properties of the parameter. This properties can be: 

• Real, Value: True or False, Default: False 

Defines, if a parameter is by Standard assumed to be real 

• Form, Value: None, Diagonal or Scalar, Default: None 
See 11221 

• LaTeX, Value: None or WT^name 

The name of the parameter used in the I^TjT;Xoutput. Standard I^T]T;Xlanguage, but for each \ you 
have to use \\. 

• Dependence, Value: None or Function 

The parameter is always replaced by this definition, see 14.2.31 

• DependenceOptional, Value: None or Function 

It can be chosen, if the parameter is replaced by the made definition. see 14.2.31 

• DependenceNum, Value: None or Function 

This definition is used in numerical calculations, see 14.2.31 

• MatrixProduct, Value: None or List of two matrices 

The parameter is defined as a product of two matrices, see 14.2.31 

• LesHouches, Value: None or Position in Les Houches file 

The numerical value of the parameter can be given by a Les Houches file, see 14.2.31 

• Value, Value: None or Number 

The numerical value of the parameter can be defined here, see 14.2.31 



29 



The Model Files 



4.2.2 Simplifying Assumptions 

SARAH normally handlcs parameters in the most general way. This means parameters normally assumed 
to be complex and all matrices can have off diagonal values. This can be changed by certain statements 
in the parameter file: It is possible to define a parameter as real by setting 

Real -> True 

The gauge couplings are by default assumed to be real. 

The degrees of freedom for a parameter, which is in principle a tensor, can be reduced by using the Form 
statement witli the foUowing options 

• Diagonal: Only diagonal entries are assumed to be unequal from zero 

• Scalar The tensor is replaced by a scalar 

Examples CP and Flavor conserving Yukawa matrices are defined by 

{Yu, { Real -> True, 

Form -> Diagonal}} 

Another example are the gaugino masses. Normally, all generations of gauginos have the same mass, so 
the parameter can expressed by a scalar instead of a matrix: 

{MassG, -C Real -> True, 

Form -> Scalar}} 

4.2.3 Defining Values and Dependences for Pcirameters 
4.2.3.1 Dependences 

There are different possibilities to define dependences between parameters by using the Dependence 
statements. The difference between the three statements is the time at which the relations are used 

• Dependence: The relations are always used. The correspondingsubstitutions aresavedin subAlways. 

• DependenceOptional: This relations are only used, if the options UseDependence is set to True, 
e. g. by calculating vertices. The substitutions are saved in subDependences. 

• DependenceNum: This dependences are only used, if a numerically value for the parameter is calcu- 
lated. The substitutions are saved in subNum. 

The indices of vectors or tensors are implicitly assumed to be indexl, index2, This in combination 

with sum[index, start, final] can be used in the foUowing way 

{X, {Dependences -> sum[nl,l,3] suin[n2,l,3] A[indexl,nl] Y[nl,n2] B [n2 , index2] }} 

is interpreted as 

3 3 

-^iii2 ^ ^ ^ ^ ^ ^iirai^in2-^n2i2 (4T) 
ni=l 712 = 1 

Parameters can also defined as matrix product of other matrices by using MatrixProduct. The argument 
must be a list consisting of two matrices of same dimension: 

{X, {MatrixProduct-> {A, B} }} 

With this definition, every appearance of a matrix product of A'^B is replaced by X and B^A by X^ . 
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Examples 

• A dependence always might be used is the parametrization of a mixing matrix by a inixing angle: 
The niixing of the charged Higgs (ZP) iii the MSSM is parametrized by a inixing angle /3. This can 
be done in SARAH with 

{ZP, i Dependence -> {{-Cos [\ [Beta] ], -Sin[\ [Beta] ]} , 

{-Sin [\ [Beta] ] , Cos [\ [Beta] ] » 

• An exaniple for an optional statement might the relation between the gauge couphngs gi and 172 
and the electric charge e. This relation is defined by 

{gl, {DependenceOptional -> e/Cos [ThetaW] }} 

Now the result for vertices can be expressed in terms of the electric charge by 
Vertex[List of Particles, UseDependences -> True] 

• A relation which might only be used for numerical calculations is the relation between the Weinberg 
angle and the gauge couplings: 

{ThetaW, {DependenceNum -> ArcCos [g2/Sqrt [gl~2+g2~2] ] > 

• The CKM matrix is defined as the product of two rotation matrices: 
{CKM, {MatrixProduct -> {Vd,Vu} }} 

4.2.3.2 Numerical Values 

If the considered parameter does not depend on other parameters, there are two ways to assign a numerical 
value to this parameter: 

• Value: Adds directly a numerical value to the parameter definition 

• LesHouches: Defines the position of the numerical value for the parameter in a Les Houches spec- 
trum file (see also |18|). 

Example The numerical values of the CKM matrix are known, and can be defined as 



> {{{1,1}, 


0. 


,97383}, 


{{1,2}, 


0. 


,2272}, 


{{1,3}, 


0. 


,00396}, 


{{2,1}, 


0. 


,2271}, 


{{2,2}, 


0. 


,97296}, 


{{2,3}, 


0. 


,04221}, 


{{3,1}, 


0. 


,00814}, 


{{3,2}, 


0. 


,04161}, 


{{3,3}, 


0, 


,999100}} 



}} 
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The values for the soft breaking parameters are normally calculated by using a spectrum calculator like 
SPheno [20] and written in a Les Houches file. To read this file with SARAH, it is possible to define the 
position for the values in the spectrum file with 

{md, { LesHouches -> iiil,!}, ■[MS0FT,47», 

{{2,2}, {MS0FT,48}}, 
{{3,3}, {MS0FT,49}}} 

}} 

4.3 Particles File 

The particle file contains Information about the fields of the model, which are needed for the out- and 
input. 

• RParity: The R-Parity of a particle: +1 or -1. If not defined, +1 is used. 

• PDG: The PDG number: Needed, if the written model files should be readable by event generators 
or if masses are given by a Les Houches file. If not defined, 99 will be used. 

• Width; The width of the particle. If not defined, O is used. 

• Mass: The options of defining the mass of a particle are: 

— Numerical Value: A numerical value for the mass of the particle is given. 

— Automatic: SARAH derives the tree level expression for the mass from the Lagrangian and 
calculates the value by using the values of the other parameters. 

— LesHouches : SARAH reads the mass from a Les Houches file. 

• Latex; The name of the particle in I^T[T;Xfilcs in Standard I^T[T;Xlanguage. But for each \ one has 
to use \\. If not defined, the Mathematica InputForm of the particle name is used. 

• FeynArtsNr: The number of the particle in a FeynArts model file. If not defined, the number will 
be generated automatically 

• Output: A short form the particle name consisting of two letter with no no-standard signs. Needed, 
to make sure that all programs outside SARAH can read the name correctly. If not defined, the 
Mathematica InputForm is used. 

This Information must be given for all eigenstates in arrays named 

ParticleDef inition ["Name of Eigenstates"] 

Only for the I^T[;]Xoutput also the names of the Weyl spinors and intermediate states (like scalar and 
pseudoscalar components of Higgs) should be given to improve the layout of the produced pdf file. 

Example For the eigenstates after EWSB a entry might look like 
ParticleDef initions [EWSB} = { 

{Sd , { RParity -> -1, 

PDG -> {1000001,2000001,1000003,2000003,1000005,2000005}, 
Width -> O, 
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Mass -> Automatic, 
FeynArtsNr -> 14, 
LaTeX -> "\\tilde{d}" , 
OutputName -> "dm" }}, 
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Chapter 5 

Defining Models 



5.1 Particle Content 
5.1.1 Vector Superfields 

First you have to define the gauge sector of your model. This is done in the following way: Every entry 
of the array Gauge stands for vector superfield and contains the foUowing information: 

Gauge [ [i] ] ={Superf ield Name, Dimension, Nsune of Gauge Group, Coupling, Expand}; 

The different parts have the fohowing meaning: 

• Superfield nsune: 

This is the name for the complete vector superfield containing vector bosons and gauginos. This is 
also the basis of their name as explained in 13.11 

• Dimension: 

This defines the dimension of the SU{N) gauge group: U[l] for an abelian gauge group, or SU [N] 
with integer N for a non-abelian gauge group. 

• Ncune of Gauge Group: 

This is the name of the gauge group, e. g. hypercharge, color or left. This choice is import, because 
all matter particles charged under a non abelian gauge group do carry an corresponding index. The 
name of the index is the first three letter of the name plus a number. So it must be taken care that 
the first three letters of an gauge group name are not identical 

• Coupling: The name of the coupling constant, e. g. gl 

• ExpEind: Values can be True or False. If it is set to True, all sunis over the corresponding indices 
are evaluated during the calculation of the Lagrangian, see (jS.l.l.O.ip . 

5.1.1.0.1 Expand Gauge Indices Gauge indices for non-abelian gauge groups are not expanded, if 
the corresponding gauge group remains unbroken (e. g color group in MSSM). Before for broken groups, 
it is convenient to expand this indices (e. g. SU{2)l in MSSM). This can be parametrized in the following 
form: 




(5.1) 
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The sums are over the color indices a and the isospin index i. 

S ARAH is not restricted to the case, that only the indices of one gauge group are extended, but different 
or even all gauge groups can be expanded. This enables for example the possibility to study SU (2) r x 
SU{2)l. 

5.1.1.0.2 Gaugino Masses SARAH adds for every vector superfield a soft breaking gaugino mass 
Mass<>"Superf ield Name" 

For non abehan gauge groups, this mass does also carry an generation index, but in the parameter file 
(see I4.2p it can be chosen that ah generations have the same mass. 

Example: Standard model color group 
Gauge [[3]] = {G, 3, color, g3, False}; 

The consequence of this entry is 

• Gluon and Gluino are named VG respectively f G 

• The SU{3) generators, the GeU-Mann matrices, are used 

• The color index is abbreviated colX (for X = 1,2, ...) 

• The strong coupling constant is named g3 

• The sums over the color indices are not evaluated 

5.1.2 Chiral Superfields 

The matter fields are arranged in chiral superfields. 

Field[[i]] = {Components , Generations, Superfield Name, Transf ormation 1, Trainsf ormation 

• Components: The basis of the name for the components. Two cases are possible: 

— The field transforms only trivially under the gauge groups with expanded indices. In this case, 
the entry is one dimensional. 

— The field transforms non-trivially under gauge groups with expanded indices. In this case, the 
entry is a vector or higher dimensional tensor fitting to the dimension of the field. 

• Generations: The number of generations 

• Superfield Nsune: The name for the complete superfield 

• Transf ormation 1, Transf ormation 2,..: Transformation under the different gauge groups de- 
fined before. For U{1) this is the charge, for non abelian gauge groups the dimensions is given as 
integer respectively negative integer (any irreducible representation is possible, see|B]). There must 
be so many charges as defined gauge groups! 

Dimension D not unique! The dimension D of an irreducible representation is not 
necessesarily unique. Therefore, to make sure, SARAH uses the demanded represen- 
tation, also the corresponding Dynkin labels have to be added! 
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Examples 

a) Fields with expanded indices 

As exainple the definition of the left handed quark superfield, if you assume only the Standard 
model gauge groups: 

Field[[l]] = {{uL, dL}, 3, q, 1/6, 2, 3>; 

The consequence of this definition is 

• Left handed up-squarks and quarks are called SuL / FuL 

• Left handed down-squarks and quarks are called SdL / FdL 

• There are three generations 

• The superfield is named q 

• The soft breaking mass is named mq2 

• The hypercharge is i 

• The superfield transforms as 2 under SU{2) 

• The superfield transforms as 3 under SU{3) 

b) Fields with no expanded indices 

As example the right handed down quark superfield is defined in the MSSM as 

Field[[3]] = {{conj[dR]}, 3, d, 1/3, 1, -3}; 

This means the following: 

• The right handed squarks and quarks are called SdR and FdR 

• There are three generations 

• The Superfield name is d 

• The soft breaking mass is named md2 

• The hypercharge is i 

• It does not transform under SU{2) 

• It does transform as 3 under S't/(3) 

c) Speficiation of representation. Since, the 10 under SU (5) is not unique, you can define it, using the 
corresponding Dynkin labels. 

Field[[l]] = {Ten, 1, t, {10 , {O , 1 , O , 0}} , . . . } ; 
or 

Field[[l]] = {Ten, 1, t, {10 , {O , O , 1 , 0}} , . . . } ; 

Soft breaking mass SARAH adds automatically for all chiral superfields soft breaking squared masses 
named 

m o "Name of Superfield" <> 2 

Generally, this are quadratic, complex matrices. If only diagonal entries or no CP violation should be 
assumed, this can be done by statements in the parameter file fsee 14. 2|) . 
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5.2 Superpotential 



The definition of the superpotential is straight forward: 



SuperPotential = {-[-[Coeff icientjPareuneter, (Contraction)}-, 

{Particle 1, Particle 2, (Particle 3)} }, ...} 



Each term of the superpotential is defined due to two list: 

• The first list is two dimensional and defines a numerical coefiicient (Coef f icient) and the nama 
of the parameter (Parameter) 

• The second list consists of two or three particles for the bi- respectively trilinear terms 

• The contraction of the indices can be given optionally. 

Contraction of indices The indices of the involved particles are automatically contracted by SARAH 
using Kronecker deltas and epsilon tensors. Sometimes, there are more possibilities to contract all indices. 
Therefore it is also possible to iix the contraction of each term. 
The contraction used by SARAH can be seen by using 

ShowSuperpotentialContractions ; 

Properties of Couplings and soft breaking terms If the particles involved in the different inter- 
actions have more than one generation, the couplings are in general complex tensors carrying up to three 
generation indices. Assumptions like diagonality or no CP value can add by using the parameter file (see 



The corresponding soft breaking term to each superfield coupling is automatically added to the La- 
grangian. The soft breaking couplings are named as 

A["Name of Coupling"] 

and carry the same indices as the superpotential coupling. Simplifying assumptions can be made inde- 
pendently of the assumptions for the superpotential parameters in the parameter file. 

Example The term involving the up Yukawa coupling is 



This can be defined in SARAH with 

{{1, Yu},{q,Hu,u» 
The explicit contraction of the indices would read as 

{{1, Yu,Delta[coll,col3] Eps [lef 1 , lef 2] },{q,Hu,u}} 

There is a soft breaking term automatically added to Lagrangian, which has the same meaning as 



1121). 



Y 

' r. 




(5.2) 




(5.3) 
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5.3 Properties of different Eigenstates 

For defining the properties of the different sets of eigenstates, the DEFINITION statement is used: 
DEFINITION["Name of States"] ["Property"]= {...}; 

The possible properties are VEVs, GaugeSector, MatterSector,GaugeFinxing and Additional. 

But, boforc defining this properties, the names for all eigenstates must be fixed in the correct order. This 

is dono duo to 

NameOf States={List of Name}; 

The hst giving the names can in principle arbitrary long. Common entries are, e. g. 

NameOf States={GaugeES, EWSB> 
NameOf States=-CGaugeES, SCKM, EWSB} 

5.3.1 Decomposition of Fields 
5.3.1.1 Introduction 

The particles responsible for breaking a gauge symmetry receive a VEV. After the symmetry breaking 
this particles are parametrized by a scalar and a pseudo scalar a part and the VEV v: 



5.3.1.2 Implementation in SARAH 

This is in SARAH done by 

DEFINITION ["Name of Eigenstate"] [VEVs] = 
{Particle Name, {{VEV, Coefficient 1}, 

{Pseudoscalar , Coefficient 2}, {Scalar, Coefficient 3}}; 

• Name: The name of the particle receiving a VEV 

• VEV: Name of the VEV 

• Scalar: Name of the scalar component 

• Pseudoscalar: Name of the pseudo scalar component 

• Coefficient 1,2,3: The different (numerical) coefficients. 

All indices carried by the particle receiving the VEV are automatically added to the scalar and pseudo 
scalar part. The scalar, pseudo scalar and the VEV are handled as real parameters in SARAH. 

Example In the MSSM the Higgs gets an VEV u„ and is therefore parametrized by 



S = —={4>s^ i<ys + Vs) 



(5.4) 



= {Vu + i<Ju + (pu) 
This is done in SARAH by using 



(5.5) 



DEFINITION [EWSB] [VEVs]= 

{{SHdO, {vd, 1/Sqrt[2]>, {sigmad, \ [Imaginaryl] /Sqrt [2] } , {phid, 1/Sqrt [2] }}> ; 
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5.3.2 Mixings in the Gauge Sector 
5.3.2.1 Introduction 

After breaking a gauge symmetry, the vector bosons mix among each other as well as the gauginos do. 
Because, in general, the different generations of a vector boson or gaugino rotate to difFerent mass eigen- 
states, it is not possible to use the same parametrization as in the matter sector, shown in 15.3.31 



5.3.2.2 Implementation in SARAH 

It is possible to define the mixing in the gauge sector in the foUowing way: 
DEFINITION["Name of States"] [GaugeSector] = 

{{Vector Boson, {Generation 1, {{New la, Coeff la}, {New 2a, Coeff 2a},..}, 
{Generation 2, {{New Ib, Coeff Ib}, {New 2b, Coeff 2b},..}, 
... }, 

{Gaugino, {Generation 1, {{New la, Coeff la}, {New 2a, Coeff 2a},..}, 

{Generation 2, {{New Ib, Coeff Ib}, {New 2b, Coeff 2b},..}, 
... }}; 

Here New XY is the namc of the new mass eigenstates and Coeff la is the mixing eleinent of the new 
and the old eigenstates. SARAH interprets this definition as a suin of sums: 

N^-l / \ 
Vi^ ^ (5n, i ^ CnmVnm (5.6) 
n—l \ m— 1 / 

V is the former eigenstate with generation index i running from 1 to N"^ — 1, c„m are the coefficients for 
the new eigenstates Vnm- 

Example We consider the electroweak symmetry breaking of the SM. The three W-Boson Wi mixes 
to charged eigenstates and W~ = W'^ and one part of the photon 7 and the Z Boson Z: 



Wi 



L[w-+W-*) (5.7) 



n/2 



W2 = {w-* - W-^ (5.8) 

Wy. = sinQwl + cofiQwZ (5.9) 
This is done by 

{VWB, {l,{VWm,l/Sqrt[2]},{conj [VWm] , 1/Sqrt [2] }} , 

{2 , {VWm, -\ [Imaginaryl] /Sqrt [2] } , {conj [VWm] , \ [Imaginaryl] /Sqrt [2] }} , 
{3,{VP, Sin[ThetaW]},{VZ, Cos [ThetaW] }}} , 

in SARAH. 



In the definition of electroweak symmetry breaking (VWp) can not be used! 
This would introduce a new complex field not related to W~ (VWm) and therefore 
change the degrees of freedom of the theory! Of course, the winos and W~ 
are not related by complex conjugation, so two new fields f Wm and f Wp are used 
after mixing. 
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S ARAH always handles the vector bosons for unbroken gauge theories as real particles. If the definition 
of the new mass eigenstates is free of coinplex parameters (hke in the case of photon and Z) SARAH 
handles also this new eigenstates as real parameters. 

5.3.3 Particle Mixings 

Symmetry breaking or also bilinear terms in the superpotential lead to a rotation of the former (gauge) 
eigenstates to new mass eigenstates. How to define this rotation in the model file depends on the fact, if 
the mass matrix is symmetric or not. 

5.3.3.1 Introduction 

Properties of symmetric Mass Matrices In the symmetric case the Lagrangian has the form 

^Mass = (I^'^MC^. (5.10) 

The matrix M can be diagonalized by a matrix U: 

Mma = U-HlU. (5.11) 

The eigenvalues of the mass matrix M are the masses of the new mass eigenstates -0, which are related 
to the former eigenstates by 

V'' = V<j)j (5.12) 

Properties of Non-Symmetric Mass Matrices If the mass term in the Lagrangian is built by two 
vectors (f)i,(f)2, i.e. 

L = (l)lM<j)2 (5.13) 
two niixing matrices are needed to diagonalize the mass matrix: 

MD^a = V-^IU. (5.14) 

The two mixing matrices V and U diagonalize the matrices MM'^ and M'^M 

Af2 = V-H'IM^V, Ml^^ = U-'^M^MU (5.15) 

and connect the new eigenstates tpi and tp2 with the old ones by 

ri^V^jcl^i, 4^l^U,,<j>{ (5.16) 

5.3.3.2 Implementation in SARAH 

Both type of mixings are defined due to 

DEFINITION ["Name of States"] [MatterSector] = { ... } 

SARAH differs between symmetric and non-symmetric mixings by the form used for the definition. Note, 
that first all generations of one particle in the basis are inserted, before the next particle follows, see 
I|5.17ll5.18p . 
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5.3.3.2.1 Symmetric Mixings in SARAH The form for defining symmetric rotations is 
{{List of Old Eigenstates} , {Name of New Eigenstates, Name of Mixing Matrix}} 

• First a list of the names of old eigenstates is given. There are no indices allowed. 

• The basis for the name of the new eigenstates must be given. 

• The name of the mixing matrix, which appears in the vertices, must be given. 

Examples The mixing in the down-squark sector is given by 
{{SdL, SdR}, {Sd, ZD» 

This means, the three left handed d-squarks SdL and the three right handed d-squarks SdR are mixed to 
new eigenstates called Sd with generation index running from 1 to 6, and the mixing matrix is called MD. 
The above statement leads to the following relation between the states 

d' = E ^.^-4 (5-17) 

d'+' = E^.>3,,4 (5.18) 

In the fiavor conserving case, this matrix is of course reducible. Therefore, SARAH checks all matrices 
for reducibihty and sets the non block elements automatically to zero. 

5.3.3.2.2 Non-Symmetric Mixings in SARAH Non-symmetric rotations are defined in SARAH 
by 

{{{First Basis} , {Second Basis}} , {{First States, First Matrix} , {Second States, Second Matrix}}} 

Let us clarify this conventions again by an examples. 

Example We consider the chargino sector in the MSSM. This mixing is specified by 
{{{fWm, FHdm}, {fWp, FHup}}, {{Lm.Um}, {Lp.Up}}} 

This means that the gauge eigenstates W" (f Wm) and (FHdm) are mixed to the negative charged mass 
eigenstates A~ (Lm), while (fWp) and (FHup) form the new eigenstates (Lp)- The new and 
old eigenstates are connected by the mixing matrices M" (Um) and A/+ (Up). 

5.3.4 Gauge Fixing Terms and Ghost Interactions 
5.3.4.1 Introduction 



As explained in IA. 51 the general form a gauge fixing term in i?j-gauge is 
with some gauge fixing functions /°. 
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5.3.4.2 Implementation in SARAH 

SARAH calculates the ghost interactions from the gauge fixing terms. It is possible to defiiie the gauge 
fixing part of the Lagrangian for the different eigenstates by 

DEFINITION["Name of States"] [GaugeFixing] = 
{{Function, Prefactor}, ... }; 

Function is the / of (I5.20p . and the correspondmg faetor Prefactor. If the gauge fixhig functions involve 
derivatives of gauge bosons 

Der ["Gauge Boson"] 

has to be used. For Goldstone bosons, the internal name in SARAH is niust be used. 
The corresponding gauge fLxing parameters must be named 

RXi ["Gauge Boson"] 



Examples 

a) The gauge fixing term for the color group in gauge is: 

^GF = -:^\df.9\' (5.21) 
And is given in the model file by 

DEFINITION [GaugeES] [GaugeFixing] = 
{{Der[VG], -1/(2 RXi [VG] )},...} ; 

b) The gauge fixing term of Z-Boson after EWSB is fsee I A. 2 51) . 

^GF^-^ {d^Z^ + CzMzG"^ (5.22) 

The corresponding Goldstone boson is in SARAH the first generation of the CP-Odd Higgs (see 
lE.ll) . also Ah[{l}] . Therefore, we can write the gauge fixing term as 

DEFINITION [EWSB] [GaugeFixing] = 

{■[Der[VZ] + Mass [VZ] RXi [VZ] Ah[-[1}], - 1/(2 RXi [VZ] )},...} ; 



5.3.5 Additional couplings or Redefinition of existing Couplings 
5.3.5.1 Introduction 

Sometimes it might be necessary to define interactions in the Lagrangian, which can not be derived 
from the superpotential or the kinetic interaction. For example, this is the case in theories with bilinear 
R-Parity violation: The term 

^h.lHJ* + h.c. (5.23) 

is allowed by gauge invariance. 

On the other hand side, it might be necessary to change the properties of some vertices by hand. For 
example, integrating out particles might most likely spoil SUSY and therefore supersymmetric relations 
are not longer valid: the Standard model Higgs self couplings are a free parameter, while they are fixed 
in SUSY by the gauge and Yukawa couplings. 
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5.3.5.2 Implementation in SARAH 

Both, defining new interactions and changing existing ones, is done in SARAH with one statement. For 
each set of eigenstates of the model, new/changed terms can be defined separately by declaring it in the 
model file via 

DEFINITION["Name Of Eigenstates"] [Additonal] = 

{ {-CPield 1, Field2, ...}, {Coef f icient , Coupling} , Options} , ... }; 

The complete object is written directly to the Lagrangian, so it niust has mass dimension 4. First, a list 
of the involved fields is given. The second part of the definition is a numerical coefhcient and the name of 
the coupling. The couplings are handled in a similar way as the couplings of the superpotential: Tensor 
indices are added automatically, and they are assumed to be complex. Further assumptions about the 
coupling can be made in parameters .m. 

The manually defined terms will be handled like every other term of the Lagrangian, i. e. it is effected by 
rotations and replacements if the eigenstates are changed. 
The two possible options are: 

• NoHC: If this is added to a term, the hermitian conjugated of this term is not added to the Lagrangian 

• Overwrite: If this is added to a term, former couplings involving the same fields are overwritten. 
If not, the new terms are added to existing ones. 

Example 

a) Define new terms: 

The term of (|5.23p is added to the Lagrangian of the gauge eigenstates by 

DEFINITION [GaugeES] [Additional] = 

{{{SvL, conj[SHuO]>, {1, mHul2» , {-[SeL , conj[SHup]}, {1, mHul2}»; 

b) Adding additional terms to existing couplings: 
With 

DEFINITION [EWSB] [Additional] = {{{hh,hh,hh,hh}, {l.kappa}}}; 

the Higgs self couplings receive a additional contribution: 

r„4 {a,g} + 6, y/) ^ r^4 (a.g^ + b,Y^ + (5.24) 

c) Overwriting existing terms: 

To overwrite former expressions for the Higgs self interactions, use 

DEFINITION [EWSB] [Additional] = 

{{{hh,hh,hh,hh}, {l.kappa}, {Overwrite», {{hh.hh, Ah, Ah}, {1, lambda} , {Overwrite}» ; 

This has the following effect: 

r„4 {a,g1 + hY^) ^ r^4 «: (5.25) 

d) For interactions involving derivatives, use Der: 
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DEFINITION [EWSB] [Additional] = {{{Der [VG] ,conj [Se] ,Se}, {1 ,kappa}}} ; 

e) Fixing the index structure: 

It is also possible to add Kronecker deltas to the coefficient: 

{{Sd.conj [Sd] ,Se,conj [Se]}, {Delta [coli, col2] , \[Gainma]}> 

is interpreted as 



5.4 Effective or non-SUSY theories 

It is easy in SARAH to integrate particles out of the spectrum to gct an efFective thcory, or just to dclote 
this particle to get a non supersymmetric hmit of the theory. We start with the definition of effective 
theories. 

5.4.1 Integrate out Particles 
5.4.1.1 General 

If in a theory very heavy particles $ exist, this are not physical degrees of freedom if the energy scale 
is below the mass of the particle. In this case, an effective theory for the lighter fields (f) is derived by 
integrating out the heavy states 



This procedure will lead to higher dimensional opcrators like the four fermion interaction in Fermis fa- 
mous theory. In SUSY theories also other Operators like six scalar interactions are possible depending 
on the integrated out fields. 



5.4.1.2 Implementation in SARAH 

Integrating out particles in SARAH is easy: The heavy particles, which should be integrated by, are 
assigned by 

IntegrateOut[[l]] = {Particle Name} 

in the mode file. Here Particle Name is the component name of your particle, e. g SdR for all right 
handed d-squarks. If only specific generations of one field should be integrated out, the assignment is 

IntegrateOut[[l]] = {Particle Name[{first,last}]} 

Here f irst is the first generation, which should be integrated out, and last the last one. 



(5.26) 




(5.27) 



The array IntegrateOut can consist of different eigenstates of the model. All 
particles will be integrated out at the first time they do appear. 



44 



5.5. DEFINE DIRAC SPINORS 



Example 

a) In focus point regions of the mSugra parameter space the Gluino is very heavy. To get an efFective 
theory by integrating out the gluino, use 

IntegrateOut [[1]] = f G; 

b) RGE running of soft breaking masses leads to the effect that the first two generations of sleptons are 
heavier than the third oiie. Therefore, building an efFective theory with only the third generation 
of squarks is done by 

IntegrateOut [[1]] = SdR[-[l,2}]; 

IntegrateOut [[2]] = SuR[-[l,2}]; 

IntegrateOut [[3]] = SdL[-[l,2}]; 

IntegrateOut [[4]] = SuL[-[l,2}]; 

5.4.2 Delete Particles 

Deleting particles can be done in the same way as integrating them out. The difference is that there are 
no efFective operators calculated. Deleting is therefore of course faster and should be used, if the higher 
dimensional operators are not needed. 
Deleting particles is done by 

DeleteParticle={List Of Particles}; 

The properties and usage of DeleteParticles is the same as of IntegrateOut. 

What might be done with 'DeleteParticles' By deleting particles it is easy get a non supersym- 
metric limit of a model: So it is possible, to get the Standard model by deleting all SUSY particles and the 
non-standard Higgs. Since SARAH supports are second SU{3) it might also be possible to get technicolor 
(-like) theories with this procedure. 

Of course, one must be careful: Some remaining interactions are still based on SUSY relation, which are 
not longer valid for non-supersymmetric theories. For example, the four point self-interaction of the light 
Higgs is in SUSY defined by the gauge interactions, but a free parameter in the Standard model. To 
delete such relations the option for redefining couplings can be used, see 15.3.5.21 

5.5 Define Dirac Spinors 

Event generators and programs for calculating Feynman Diagrams are normally written for Dirac Spinors, 
but SARAH does all internal calculations for Weyl Spinors. Therefore it is necessary to define, how the 
Weyl Spinors combine to Dirac Spinor by 

Dirac = {Weyl 1, Weyl 2} 

Dirac is the new name of the Dirac spinor, while Weyl 1 and Weyl 2 the names of the Weyl spinors 
building the left respectively right component of the Dirac spinor. 

Example The electron Fe is built from the components FeL and FeR by 
Fe = {FeL, FeR} 

while the neutralinos Chi, Majorana particles, consists only of the mass eigenstates LO 
Chi = {LO, conj [LO] } 
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How to generate output files by using the command ModelOutput was already explained in section[X71 
Here we want to give some further information what the different routines do actually. 

6.1 LaTeX 

6.1.1 Writing a LaTeX File 

It is possible to write a MJfpCfile with all information about the model you are considering by using 
ModelOutput [Eigenstates , WriteTeX->True] ; 

If you have aheady calculated aU vertices by using ModelOutput, but not generated a M^rjKfile, you can 
directly generate the TeX output with 

MakeTeK [] ; 

This uses the existing results and saves therefore time. 

There are different tex-files produced containing the following information: 

• List of the fields 

• Mass Matrices, if there is any mixing 

• Renormahzation Group Equations 

• Tadpole equations, if there is a broken gauge symmetry 

• One foop self energies and tadpole equations 

• All interactions 

SARAH will translate all internal names to the defined I^TfTpCoutput by using string replacements with 
regular expressions. The reason for this is that we want to be compatible also with Mathematica 5.2 and 
there are some difficulties in defining a TeKForm for symbols for which already a DisplayForm is defined. 
Also we were not completely satisfied with the Standard TeX output of Mathematica. Therefore we have 
not just used the Mathematica internal command TeXForm to generate the M^rpC-code. 
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6.1.2 Options 

There are several options for the M^rjKoutput, which can be set with 
SetOptions [WriteTeX, Options] ; 

permanently for the Mathematica session, or by appending it to the write command 
MakeTeK [Options] ; 
The options are 

• FeynmanDiagrams, Values: True or False, Default: True 

If True, there will a Feynman Diagram generated for each vertex. This is done by using the I^T[5X- 
Package FeynMP from Thorsten Ohl [21], which can be found in the directory LaTeX-Packages it is 
not installed yet. 

• Ef f ectiveOperators, Values: True or False, Default: True 

If analyzing an effective theory, there will also higher dimensional operators exist. With this option 
it can be chosen, if also the vertices for the higher dimensional operators are included in the 
I^TgKfile. By default, this are only the vertices involving up to four particles. For switching on six 
particle interactions,SixParticleInteractions is used. 

• SixParticleInteractions, Values: True or False, Default: False 

If a scalar or a vector boson is integrated out, also six particle interactions of scalars and vector 
bosons may be generated. This vertices are often lengthy and the calculation of this vertices will 
take seriously longer then for three and four particle interactions. Therefore this statement is 
included to give an explicit choice, if this vertices are calculated and if they should be inserted in 
the tex file. 

6.1.3 Making the pdf File 

The tex files are saved in the directory 

. ./Output/"Name of ModelV'Name of Eigenstates"/TeX 

and the main file is "Ncune of Model" <> "Name of Eigenstates" o.tex. Ali other files are included 
in this file by using the input-command of WT^^^. If Diagrauns is set to True, the foUowing steps must 
be done for generating an pdf document including the diagrams: 

a) First conipile the tex file, e. g. pdflatex model.tex 

b) Go to the directory Diagrams and compile every .mp file with mpost. This is done under Linux and 
under Windows with 

mpost FeynmcLnDiaX.mp 

It is also possible to apply the mpost command on all .mp-files at once by using 
find . -nsune "*.mp" -exec mpost {} \; 

c) After generating all diagrams, go back and compile the .tex-file again by using pdflatex. 
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To simplify your life a little bit SARAH will write a shell script in the tex-output directory, which does 
exactly this three steps. It can be started uiider Linux with 

. /MakePDF . sh 

or under Windows with 
MakePDF . bat 

It is possible that the script must be first declared is executable in Linux by 
chmod 755 MakePDF. sh 

6.1.4 m^Output for new Models 

If you have defined your own model, it is suggested to include the M^rjK-code for all names of the new 
fields and paraineters in parameters .m and particles .m. For more information, see 14. 31 and 14.21 

6.2 FeynArts 

6.2.1 Generate a FeynArts Output 

A model file for FeynArts is created by 
ModelOutput [Eigenstates , WriteFeynArts->True] 

or, if ModelOutput was used before and the vertices are already calculated, 
MakeFeynArts 

The following things are done: 

• A list of all particles in the models is generated. The particles in FeynArts are named 

— S [X] : For scalars, with some integer X 

— F [X] : For fermions, with some integer X 

— V [X] : For vector bosons, with some integer X 

— U[X] : For ghosts, with some integer X 

X can be defined in the particle definition file of SARAH or is chosen automatically. FeynArts also 
supports labels for particles which are easier to read for humans eyes by using a TeX-like output. 
The label for each particle is therefore generated from the defined M^TjKname in SARAH. 

• A list with all appearing indices is written. This are indices for generation and non abelian charges. 
In contrast to the originally model file for MSSM of FeynArts, every particle gets an independent 
generation index. 

• The list with interactions is written. If the theory contains more non abelian, unbroken gauge 
groups, the generators of this gauge groups will appear in the vertices, of course. By default, the 
generators of SU (3) are associated to SUNT and so automatically simplified by using FormCalc |22| . 
The algebra for generators of another dimension must be done separately. 
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6.2.2 Dependences, Numerical Values and special Abbreviations for Form- 
Calc 

It is a second file generated by SARAH together with the model file for FeynArts: Substitutions-<> 
"Najne of Eigenstates" O. m . This file contains additional information whicli might be useful for 
calculating diagrams: 

• Replacement rules with tlie defined dependences in SARAH, see 14.2.31 

— Dependences 

— DependencesOptional 

— DependencesNum 

• The definitions of the masses: MassesSARAH 

• The numerical values for the parameters in SARAH: Numerical Values 

• Special abbreviations for FormCalc as they are also defined for the MSSM and SM in the FormCalc 
package: 

— A complex conjugation is replaced by C: Conjugate[X] -> XC 

— A square is replaced by 2: X~2 -> X2 

— The names of soft breaking couplings are merged: A[X] -> AX 

It is suggested to use this definitions to speed up the calculations with FormCalc. 

6.2.3 Higher Dimensional Operators 

The Standard version of FeynArts does not support higher dimensional operators, because the corre- 
sponding Lorentz structures are not given in Lorentz . gen. Therefore, such operators are not included 
in the model file, even if an effective theory is considered. 

6.3 CalcHep/CompHep 

6.3.1 Generate Output for CalcHep and CompHep 

ModelOutput ["Name of Eigenstates", WriteCHep->True] 

In this case the Standard values for the options are used. If the vertices have been calculated already by 
ModelOutput, 

MakeCHep [options] 

can be used. The possible options are 

• FeynmanGauge, values: True or False, Default: True 

By setting to True, the interactions of the Goldstone bosons are written in the interaction file. 

• CPViolation, values True or False, Default: False 

By setting to True, the possibility of CP violation is included in the model files, see section [6.3.51 
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• ModelMr, values Integer, Default: 1 

The number added to file iiames, see next section. 

• CompHep, values: True or False, Default: False 

By setting to True the model files are written in the CompHep format. 

6.3.2 Model Files 

The CalcHep/CompHep output function of S ARAH generates the following four files 

• prtclX.mdl: Contaimng all particles 

• lgrngX.mdl: Contaimng the interactions 

• varsX.mdl: Containing the numerical values of the variables 

• funcX.mdl: Containing dependences between the parameters 
X is a number, which can be chosen by the option ModelNr. 

6.3.2.1 Particles 

First, there are stringent constraints on the naming of particles in CalcHep: Only names up to four 
letters are allowed and also a index structure is not supported. Therefore, it is not possible to use the 
SARAH internal definitions of particles. Instead the names used for the model files are based on the 
defined OutputName of each particle in the following way 

a) The basis of each name is the entry in OutputNeune in the particle file, see 14.31 

b) If the considered particle is not self-conjugated, for the anti particle the first letter is changed from 
upper to lower case or vice versa. 

c) If there are more generations for one particle, the number of generations is appended at the end of 
the name 

d) If the defined R-parity is -1, a ~ is added to the beginning of the name to assign SUSY particles. 
So it is possible to use the model files in MicrOmegas without the need of an additional list of all 
SUSY particles 

This is the Standard procedure for all vector bosons, fermions and most scalars. Ghosts, Goldstone bosons 
and auxiliary fields handled in a different way: 

There are three different kind of ghosts. These are not written in the particle file, but appear in the 
Lagrangian file: 

• Fadeev Poppov Ghost: This are the well known Ghost derived from the gauge transformations of 
the gauge fixing term. The name in the model file is 

"Name of Vector Boson".C 

• Goldstone Ghosts: This are just the Goldstone bosons 'eaten up' by the gauge particles. Their 
name is 

"Ncune of Vector Boson".f 
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• Tensor Ghosts: Is needed to express the four gluon interaction. The name is 
"Name of Gluon" .t 

SARAH derives the name of Goldstone and Faddeed-Poppov ghost automatically from the underlying 
vector boson, but the tensor ghost and its one interaction with two gluons is hard coded. 

The last kind of fields known by CalcHep/CompHep are auxihary fields. Their purpose is explained in 
the next section, but their names are as foUows 

~0X 

where X is a number. The antiparticle, if it is not the same, is counted as X+l. 



6.3.3 Auxiliary fields in CalcHep/CompHep 

We mentioned in the last section, that CalcHep and CompHep needs special auxihary fields. The reason 
for this is that the color structure is hard coded, and therefore interactions of four colored particles or 
two colored and two gluons sufFer from an ambiguity. Therefore, this interactions are split in two three 
particle interactions by inserting auxiliary fields. 

SARAH does a similar splitting for all interactions between for complex scalars by inserting auxiliary 
fields when calculating the F- and D-Terms. Also the interactions between two squarks and two gluons 
are split in two three particle interactions. 

6.3.4 Vertex Functions 

All interactions are parametrized by a variable. The value of this variable is then defined by using the 
results of SARAH for the corresponding vertices, but the following renaming had to be done: 

• Tensor indices are just added to the name, therefore all sums in the vertices had do be evaluated: 
sum[il,l,3,MD[l,il]] -> MDll + MD12 + MD13 

Some values, which are known to be zero like in the flavor conserving case, are set to zero at this 
point. 

• Variables names, which are longer then six letters, are cut 

• All parameters are assumed to be real, so complex conjugation is remove (see 16.3.5)) 

• The generators and structure constants of strong interaction are removed, because they are defined 
implicitly in CalcHep 

Since there is no index structure allowed, all interactions must be written for all possible combinations 
of generation indices. This can lead for some models to a very long time for writing the model file, e. g. 
in the fjLvSSM with 8 charged Higgs: All 8^ combinations of the self interactions must be written. 
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6.3.5 CP Violation 

CalcHep/CompHep can't handle coniplex values in the function or vars file, but only in the Lagrange 
file. Therefore all variables are by default assumed to be real, when SARAH writes the model file with 
default options. This can be changed by setting CPViolation to True. In this case, SARAH splits all 
parameters, which are not explicitly defined as real, in real and imaginary part: 

X -> RX + i*IX 

The real and imaginary part for every interaction is calculated using that splitting, and both parts are 
written separately in the Lagrange file: 

vOOOl -> RvOOOl + i*Iv0001 

6.3.6 Vciriables and Dependences 

At the end, for calculating processes, CalcHep/CompHep needs to know the numerical values for the free 
parameters and the dependences between the parameters. The numerical values are given in varsl .mdl. 
SARAH writes in this file the output name of every needed variable. Since CalcHep/CompHep does only 
support variable names with a Icngth of maximal 6 letters, SARAH cuts the name of all variables, which 
are longer then that. If a numerical value for the variable is availablc, c. g. if it is defined in the parameter 
file or it was read in from Les Houches file, it will also be added in the vars file, otherwise NaN appears. 
The dependences defined in SARAH are written as functions in funcl.mdl. 

6.3.7 What can be a problem... 

We have made the foUowing experiences by testing model files with CalcHep/CompHep: 

• A PDG number of O is not allowed for other particles than auxiliary particles. 

• In the vars file is no discrimination between small and capital letters. This must be taken into 
account by naming the mixing matrices and couplings in SARAH. 

• Higher dimensional operators are not supported 

• The color structurc is hard coded and indices are not supported in CalcHep/CompHep. So it is 
hard to implenient models with other ubroken non abelian gauge groups than the color group. 

• As already mentioned: Writing the output for models with particles appearing in a large number of 
generations (let's say, more than 6) and non reducible mixing matrices like in the flavor conserving 
MSSM, last very long, because all possible combinations of indices have to be written separately. 

6.4 Automatic Output 

It is possible, to generate the desired output automatically. To do that the the array MakeOutput in the 
model file is used 

MakeOutput = {{ Eigenstates 1 , List 1}, i Eigenstates 2, List 2}, ... }; 

Eigenstates fixes the eigenstates for which the output defined in the list should be generated. The 
entries in List X can be the foUowing: 

• TeX 
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• FeynArts 

• CHep 

Example To get a WT^- and FeynArts file for the eigenstates after electroweak symmetry breaking, 
MakeOutput = {{EWSB, -[TeX, FeynArts}}}; 
is used. 



53 



Chapter 7 

Numerical Values 



SARAH offers you some routiiie for working with numerical values of parameters and calculating mixing 
matrices and masses. 

7.1 Adding Numerical Values 

There are different ways to define numerical values for a parameter. It can just be added in the parameter 
file as described in 14.21 or the values from a Les Houches spectrum file can be used. To read in a Les 
Houches file automatically by evaluating the command Start, add 

SpectrumFile="Name of File" 

to the model file. The spectrum file must be in the same directory as the model file. 

It is also possible to read a spectrum file afterward the automatized calculation of SARAH is finished by 
ReadSpectrum ["Spectrum File"] 

An additional possibility to add a numerical value during the work with SARAH is 
SetParameterValue [parameter , value] ; 

parameter is the name of the parameter and value the numerical Value. Of course, it is also possible to 
remove easily a value. For this purpose the command 

DeleteParameterValue [parameter] ; 

is used. 

7.2 Calculate Mixing Matrices 

After numerical values for all free parameters are defined, the mass eigenstates and entries of mixing 
matrices can be calculated by 

CalcMatrices ; 

The eigenvalues and mixing matrices of each mass matrix are saved in two additional variables: 
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7.3. CALCULATE NUMERICAL VALUES 



• Eigenvalues: The eigenvalues of the mass matrix are saved in 
Mass o "Name of Eigenstates" 

• The numerical value of the mixing matrix are saved in 
"Name of Mixing Matrix" o Num 

Negative Squared Masses It could happeii, that for some parameters eigenvalues of the mass matrices 
for scalars are getting negative. If this happens, there will appear a warning in the output, and the variable 
warning will be set to True. 

Example After using CalcMatrices, the numerical value for the down-squark mixing matrix are saved 
in MDNum and the corresponding squark masses are saved in MassSd. 

7.3 Calculate Numerical Values 

To get the numerical value for a term 
Numerical Value [x] ; 

is used. X can be e. g. the entry of a mass matrix or a vertex 

Examples The numerical value for the Higgs mass matrix is calculated with 
NumericalValue [MassMatricesFullEWSB [ [4] ] ] ; 

The numerical value for the interaction between the photon and down squarks is calculated with 
NumericalValue [Vertex [{VP , Sd, conj [Sd] >] [ [2 , 1] ] ] ; 



55 



Appendix A 

Calculating the Lagrangian 



In principle, it is straightforward to get from the superpotential and gauge sector to the Lagrangian. We 
use here the method described in [23J. In brackets is always the name of the variable in SARAH given, 
which corresponds to the particular part of the Lagrangian. 

A.l Interactions of Chiral Superfields 

If we call the superpotential for a given theory W and use (p for the scalar and ip for the fermionic 
component of a chiral supermultiplet, then the interaction parts of the Lagrangian are given by: 

• Interactions from the Superpotential (Superpotential): 



-W'^i^,i>j + h.c. 



(A.l) 



• F-Terms (FTerms): 



^F"F, + h.c. 



(A.2) 



with 



54>, 



W 



(A.3) 



'J 



and 




(A.4) 



A.2 Interactions of Vector Superfields 



The most general Lagrangian only involving gauge particles (A is the spin- 
of a vector supermultiplet) is 



i and A'^ the spin-1 component 




a 



(A.5) 
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A.3. SUPERSYMMETRIC GAUGE INTERACTIONS 



with the field strength 

f;, = d^Al - d.A; + gr'"'AlAl, (A.6) 
and the covariant derivative 

D^A" = 9^A« + 5r'''=4A^ (A.7) 
This leads to the foUowing interactions: 

• Self-interaction of three gauge bosons (GaugeTri) 

= -\{d^A-^ - d,A;)gr'^A^'''A-''^ (A.8) 

• Self-interaction of four gauge bosons (GaugeQuad) 

£'AV = -\9\fabcAlAl){r''-A^'^-A'''<') (A.9) 

• Interactions between vector bosons and gauginos (BosonGaugino) 

L^y = -iAt°CT''£>^A" (A.IO) 

A. 3 Supersymmetric Gauge Interactions 

The parts of the Lagraugian with both chiral and vector superfields are 

• Kinetic Parts for Scalars (KinScalarSing+ KinScalarDoub) 

^k^n = ~D^'(t>*'D^,(t>, (A.ll) 

• Kinetic Parts for Fermions (KinFermSing+KinFerinDoub) 

^kin = -iV^'a^a^Vi (A. 12) 

• Gaugino-Fermion-Sfermion Interactions (FSGaugino): 

Lgfs = -\/2(7(0*T»A" + h.c. (A.13) 

• D-Terms (DTerms): 

Ld = ^D'^D" (A.14) 

Here are T" the generators of a gauge group. Furthermore, the covariant derivatives are 

D^<l>^ = d^cj^.-igAliT'^^), (A.15) 

D^<t>*' = d^(k*' + igAl{cj>*Ty (A.16) 

D^i}i = d^iPi - igAl{T''iP)i (A.17) 

and the D-Terms are 

D'" = -g{(j)*T''(t)) (A.18) 
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A. 4 Soft Breaking Term 

Since SUSY must be a broken theory, the following soft breaking terms can be added to the Lagrangian: 

• Soft breaking mass for every scalar (Sof tScalarMass) 

^Soft,s = m''cj)ict>* (A. 19) 

• Soft breaking mass for every gaugino (Sof tGauginoMass) 

^Soft,F = M\i\i (A.20) 

• Soft breaking terms corresponding to the terms of the superpotential 
(SoftTri+SoftBi) 

A. 5 Gauge Fixing and Ghost interactions 

Normally, the Lagrangian of a theory is invariant under gauge transformation. This invariance leads to 
difEculties in the quantization of the theory as can seen in the divergence of functional integrals. Therefore 
it is in necessary to add gauge fixing terms to break the gauge invariance. This Icads also to the necessity 
of so called Faddeev-Popov ghosts in order to give correct results for gauge invariant observables. 
The general form of the gauge fixing Lagrangian is 



£/GF 



with a function /, which can can be a function of partial derivatives of gauge bosons and Goldstone 
bosons. The corresponding ghost terms of the ghost fields fj and r] (r/, fj are not connected by complex 

conjugation), are 

'CGhost = -f]a{6n. (A.23) 
Here 6 assigns the operator for a gauge transformation. 

For a unbroken symmetry the gauge fixing terms in the so called ii^-gauge are 

^GF = -^E(^"^;)'- (A-24) 

Hcre, Vfj_ are the gauge bosons of the unbroken gauge symmetries. The theory must be independent of the 
choice of the value for the gauge fixing parameter R^, but to simplify the calculation it is also possible 
to choose a distinct value for R^. The most popular gauges are the unitary gauge — )• oo and the 
Feynman-t'Hooft-gauge R^i = 1. For broken symmetries, the gauge fixings terms are chosen in a way 
that the mixing terms between vector bosons and scalars disappear from the Lagrangian. Therefore the 
common choice for theories with the Standard model gauge sector and EWSB is 

Lgf,r, = - ^ {d'^l^f - {d'^Z^ + izMzG^f + "F^ {9^^^^ + i^^M^G+f (A.25) 

Here, and G+ are the Goldstone bosons, which make the longitudinal component of the massive vector 
bosons. 



58 



A.6. PARTS OF THE LAGRANGIAN IN SARAH 



A. 6 Parts of the Lagrangian in SARAH 

SARAH saves the different parts of the Lagrangian for the different eigenstates in the following variables: 

• LagrangianSV [Name Of States] : Parts with scalars and vector bosons (= kinetic terms for scalars) 

• LagrangianFFV[Name Of States]: Parts with fermions and vector bosons (= kinetic terms of 
scalars) 

• LagrangianSSSS [Name Of States] : Parts with only scalars (= scalar potential) 

• LagrangianFFS [Name Of States] : Parts with fermions and scalars 

• LagrangianVVV[Ncime Of States]: Parts with three vector bosons 

• LagrangianVVVV [Name Of States]: Parts with four vector bosons 

• LagrangianGGS [Name Of States] : Parts with ghosts and scalars 

• LagrangianGGV [Name Of States]: Parts with ghosts and vector bosons 

• LagrangianSSA [Name Of States] : Parts with scalars and auxiliary fields (only needed for CalcHep 
output) 

In addition, for an effective theory, there might exist 

• LagrangianSSSSSS [Name Of States]: Dimension 6 operators with only scalars 

• LagrangianSSSVVV [Name Of States] : Dimension 6 operators with scalars and fermions 

• LagrangianFFFF [Name Of States] : Dimension 6 operators with only fermions 

• LagrangianFFSS [Name Of States] : Dimension 5 operators with fermions and scalars 

• LagrangianFFVV [Name Of States] : Dimension 6 operators with fermions and vector bosons 
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Appendix B 

Group Theory 



SARAH supports not only chiral superfields in the fundamental representation, but in any irreducible 
representation of SU{N). In most cases is possible, to fix the transformation properties of the chiral 
superfield by statiug the corresponding dimension D. If the dimension is not unique, also the Dynkin 
labels are needed. 

Therefore, it is necessary to derive from this Information the corresponding generators for calculating 
kinetic terms and D-terms. Furthormore, for calculation of the RGEs the eigenvalues C2 of the quadratic 
Casimir for any irreducible representation r are needed 

T"r"(^(r) - C2{r)^{r) (B.l) 
as well as the Dynkin index / 

Tr(T'^r^)0(r) = ISab4>ir). (B.2) 
AU of that is dcrived by SARAH due to the technique of Young tablcaux. The foUowing steps are evolved: 

a) The corresponding Young tableaux fitting to the dimension D is calculated using the hook formula: 

" = n.^ (B.3) 

di is the distancc of the i. box to the loft uppcr corncr and hi is the hook of the box. 

The number of rows and columns of the Young tableaux with the correction dimension defines also 

the number of co- and contravariant indices. 

b) The vector for the highest weight A in Dynkin basis is extracted from the tableaux 

c) The fundamental weights for the givcn gaugc group are calculated 

d) The value of C2 (r) is calculated using the Weyl formula 

C2{r) = {A,A + p) (B.4) 
p is the Weyl vector. 

e) The /(r) is calculated from C2(/ ). For this step the value for the fundamental representation is 
normalized to ^. 

^(r-) = C2(r)^ (B.5) 
With -D(G) as dimension of the adjoint representation. 
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f) The number of co- and contravariant indices is cxtracted from the Young tableaux. With this 
information, the generators are written as tensor product. 

Dimension D not unique! The dimension D of an irreducible representation is not necessesarily 
unique. Especially, if you consider higher dimensional representations. Therefore, to make sure, SARAH 
uses the demanded representation, also the corresponding Dynkin labels have to be added! 
The user has access to this information: 

ChecklrrepSUN [Dim.N] 

where Dim is the dimension of the irreducible representation and M the dimension of the SU{N) gauge 
group. The result is a vector containing the foUowing information: 

• Repeating the dimension of the field 

• Number of upper indices 

• Number of lower indices 

• Value of the quadratic Casimir 

• Value of the Dynkin operator 

• Dynkin labels for the highest weight 

Examples Tho properties of a particle, transforming under the fundamental representation of SU (3) 

are calculated via 

ChecklrrepSUN [3,3] 

The output is the known result 

{3, 1, O, 4/3, 1/2, -Cl, 0}} 

The properties of a field transforming as 24 (adjoint representation) of SU{5) are calculated by 

ChecklrrepSUN [24, 5] 

The output will be 

{24, 1, 1, 5, 5, {1, O, O, 1}} 
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Appendix C 

One Loop calculations 

We list in the following the results for the generic different diagrams contributing to the one loop correction 
of the tadpoles and self energies. Different coefHcients do appear: 

• cc'. Possible charges in the loop, o. g. a factor 3 for colors. 

• Cs- Symmctry factor of ^ if thc particlcs in thc loop can not be distinguished, otherwise 1 

• Cr: Factor of 2 if the external particle is self conjugated or Majorana fermion, otherwise 1 

Cl One Loop Tadpoles 

• Fermion-Loop (FFS): 

T = 8csccmFTAo{m%) 

• Scalar-Loop (SSS): 

T = -2csccrAo{ml) 

• Vector-Loop (SVV): 

T = ecsccTAoimlr) 

C. 2 One Loop Self Energy 

C. 2.1 Fermion Field 

• Fermion-Scalar-Loop (FFS): 

E^(p2) = mFCsccCi^^)^^i■*Bo(p^m|,m|) 
= -csccCR^r^.T^/B^ip^m^ml) 

= -csccCR^TlT]:*B^{p\m%,m%) 
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C.2. ONE LOOP SELF ENERGY 



• Fermion-Vector-Loop (FFV): 

= -4csccCHmFTlr^i*Bo{p'',m''p,ml) 
E«(p2) = -csccCRrlT^£*Bi{p'',m%,ml) 
S^(p') = -csccCRT'RTl*Bi{p\m%,ml) 

C.2.2 Scalar Field 

• Fcrniioii-Loop (FFS): 

n(/) = cscccn [{rlr'^*+T],T^i*)Go{p',ml,ml) + {Tivl* +T],T^^*)Bo{p',ml,ml)) 

• Scalar-Loop (SSS): 

n(p2) = csccCRr^r^'*Bo{p'',m%,ml) 

• Scalar-Loop (4-point interaction) (SSSS): 

n(p') = -csccTAoiml) 

• Vector-Scalar-Loop (SSV): 

n(p2) = csccCRT'T^'*Fo{p'',m%,ml) 

• Vector-Loop (SVV): 

n(p') = csccCR'^T'T^'*Bo(p^,ml,ml) 

• Vector-Loop (4-point interaction) (SSVV): 

n(p2) = csccTAoim^v) 
C. 2.3 Vector Field 

We need only the diagrams involving three point interactions, because the four point interactions are 
related to them due to gauge invariance. 

• Fermion-Loop (FFV): 

n^(p2) = csccCR ((|ri|2 + \T],\^)Ho{p\mlr,m%)+4Re{TlTj,)Bo{p\m^y,m%)) 

• Scalar-Loop (SSV): 

n^(p2) = -4csccCfl|r|2B22(p',m|^,m|J 

• Vector-Loop (VW): 

n^(p2) = \Tf CSCCCR {-{4p'^+ri4^+n4jBoip'^,ml^,mlrJ - 8B22(^)^m|^,m|J) 

• Vector-Scalar-Loop (SVV): 

n^(/) = \T\''csccCRBo{p'',m^y,ml) 
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C. 3 One Loop Masses 

The one loop corrections to the masses and mass matrices can be calculated by using the results of the 
self energies. M is the corresponding mass matrix at tree level or the mass, if the fields are not rotated. 

• Majorana Fermion: 

5m = i?e(-^(5]s + Sg + S^M + IIrM + ME^ + MSl)) 

• Dirac Fermion: 

• Scalar: 

(5m2 = -i?e(n(p2)) 

• Vector Boson: 

dm^ = i?e(n(p2)) 
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Appendix D 

Gauge Anomalies 



SARAH checks if the chiral superfields of the defined models fulfill the following relations in order to 
cancel triangle gauge anomalies: 

n 

• SU{N)l. 

^Tf = (D.2) 



U{1) X SU{Nf 



2 



EE^nrr = o (D.3) 

m n 

U{1) X [/(1)2 

^^Y^i'=0 (D.4) 



• U{1) X U{l)j X U{l)k 



m n k 

The indices are for the different gauge groups are i, j, k and n, m, I for chiral superfields. 
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Appendix E 



Models 

E.l Minimal Supersymmetric Standard Model 

The name in SARAH is MSSM. 

The used conventions and parameterization are based on the ones in |24) . 
E. 1.1 Particles 

The matter fields of the MSSM are arranged in left-chiral superfields, which transform under the gauge 
groups as shown in table lE.ll The gauge fields, i. e. the gauge bosons and gauginos, are arranged in 
vector superfields transforming under the adjoint representation of the gauge groups, see table IE.21 



Name 


spin 


spin 1/2 


SU{i)c,SU{2)L.U{l)Y 


Squarks, Quarks 


Q 


{u L (11) 


[u L dL) 


(3, 2, i) 


(3 Generations) 


u 










d 


d* 




(3,1, i) 


Sleptons, Leptons 


L 


[vcl) 




(1,2,-i) 


(3 Generations) 


e 


e*R 


^R 


(1,1,1) 


Higgs, Higgsinos 


Hu 






(1,2, i) 




Hd 






(1,2,-i) 



Table E.l: Chiral Superfields of the MSSM 



E. 1.2 Superpotential 

The Matter interactions of the MSSM are a result of the following superpotential, which is a holomorphic 
function of the chiral superfields: 

W - -LY.eHd - QYddHd + QYuuHu + M^u^d- (E.l) 

Here Y,,, Y^ and Yd are the Yukawa matrices, which are 3 x 3-syinnietric matrices and often are assumed 
to be diagonal, /i is the mass term for the Higgsinos. 
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Name 


spiii 1/2 


spin 1 


SU{i)c.SU{2)L.U{l)Y 


Gluino, Gluon 


9 


9 


(8,1,0) 


Winos, W Bosons 






(1,3,0) 


Bino, B Boson 


5° 


5° 


(1,1,0) 



Table E. 2: Vector Supermultiplets of the MSSM 



E. 1.3 Soft Breaking Parameters 

SUSY is a broken symmetry, i. e. the masses of SM particles and their SUSY partners are different. This 
can be parameterized by adding so called soft breaking parameters to the Lagrangian, which introduce 
no new quadratic divergences: 

LsB = ~\'^Ma\a>^a- m^i\(t>i? 

^ A 

-{-LAeeHd - QAddHd + QAuuHu + h.c.) - {B^fiHdHu + h.c). 

They are consistmg of mass terms M a for the gauginos \a, squared mass term mf for the scalars 4>i, 
trihnear scalar coupHngs A^, often assumed to be proportional to the corresponding Yukawa coupUngs, 
and a Higgs mixing term proportional to -B^/i. 



E. 1.4 Mass Eigenstates 



The gauge eigenstates in tables lE.ll and IE.2I with the same quantum numbers can and will mix among 
each other to form the measurable mass eigenstates of the MSSM. The mass eigenstates of the MSSM are 
built up by the massive vector bosons after electroweak symmetry breaking hke in the Standard Model 
and furthermore by the neutralinos, charginos and the mixtures of sleptons, squarks and Higgs. 

E. 1.4.1 Electroweak Symmetry Breaking 

After Electroweak Symmetry Breaking (EWSB), the photon and Z-boson are mixtures of the the B- and 
Wa-boson 

Z = -sinGvi/S + cosewI^a (E.2) 
A = sinQwWz + cosQwB, (E.3) 

while the superposition of Wi and W2 form the charged eigenstates and VF^, with 

W^^-^WxTiW2. (E.4) 
v 2 

E. 1.4. 2 Neutralinos and Charginos 

In an analog way to the charged bosons the first and second Wino Wi^2 mix after EWSB to charged 
fermions 

W^ = -^WxTiW2. (E.5) 
V 2 
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Furthcrmorc. thc chargcd components of the Higgsinos, H^,H^ , and the charged m.ix to form new 

mass eigenstates, called charginos 



{Z- 



M2 

evi 



■splsw 



The four component fermions look like 



( + -- \ ^i.C _ / - -+ A 



(E.6) 



(E.7) 



with 

H+ = Zl^Kt, H^=Zl'n-, W^=iZlKt. (E.8) 

The neutral components of the Higgsinos, the Bino and the third Wino mix to four neutral mass eigen- 
states, the neutrahnos x°, which are Majorana fermions. The mass matrix in the basis (5, M^, if^, 
is 

/ Ml O 

O M2 



■'N 



2cyv 

\ 2cw 



evg etil 
— eV2 



2s 



w 



—evi 


eti2 


\ 


2cw 


2civ 


evg 


-2t)2 






2sw 







-P' 




-M 





/ 



Zn = 



V o 



/ 



and the neutrahnos read in four component notation as 



(E.9) 



Xi 



(E.IO) 



with 



(E.ll) 



E.1.4.3 Sleptons 

The three generations of the fields bl and mix to six charged eigenstates Em,i ■ ■ ■ E^fi 



2-f ali* £i 



The mixing matrix can be read of the Lagrangian as follows 

/ ^4„i ••• o 



(MDL {Ml) 



E)RR 



\ 



J 



(E.12) 
(E.13) 



(E.14) 



{MDll = 



„ o o -E' H 



4r2 



■ m 



Er 



{MDlR = -j=iV2YefX* +ViAe). 



(E.15) 
(E.16) 
(E.17) 
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E.1.4.4 Squarks 

The up-type squarks u£ and mix to tJrn,i 

ui = SljU^.^ (E.18) 

u'^ = S\j+'^"U*^,,, (E.19) 

and the down-type and to Dm,i 

di = Sj]*D,n,, (E.20) 

7/ _ n{I+3)i 



4 = S^i+'^^i)*,. (E.21) 



The corresponding niixing matrices could be read of the Lagrangian in the same way as in the case of 
the slepton mixing and are hsted e. g. in ^25] . 

E. 1.4. 5 Sneutrinos 

The three sneutrinos do mix to the mass eigenstates Vm'- 

v' = S'y'Ura.r (E.22) 

and the mass matrix is 

Ml^'M_^E + ml. (E.23) 
E.1.4.6 Higgs 

The Higgs fields will be parameterized via 



(E.24) 
(E.25) 



Here Vu and Vd are the Vacuum Expectations Values (VEVs) of the Higgs and they are connected to the 
mass of the W-boson and Z-boson 



Mw = 7^\lvl + vl (E.26) 



M z = + (E.27) 



and their ratio defines tan /3: 

tan^ = — . (E.28) 

Vu 

(pd and 4>u are the scalar components, while ad and (T„ are the pseudo scalar ones. 

Now, the scalar, pseudo scalar and charged components do mix even without EWSB and we are getting 
the following eigenstates: 
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• neutral, scalar Higgs h, H 




cos a — sm a 
sin a cos a 




(E.29) 



charged, scalar Higgs and Goldstone modes 




sin /3 — cos \ ( (t^d 
cos /3 sin /3 / \ 



with </,+ = (0-)*,0- = (^+)*. 
• CP-odd Higgs and neutral Goldstone boson G" 



(E.30) 




sin /3 — cos (i \ ( <Jd 
cos /3 sin /3 / \ (t„ 



(E.31) 



E. 1.4. 7 Gluons and Gluinos 

The gluons and gluinos don't mix. 



E. 1.4.8 Overview 

A collection of the gauge and corresponding mass eigenstates is shown in table IE.3I Often the mixture 
of the first and second generation of squarks and sleptons is negligible because of the small Yukawa 
couplings, i. e. the mass eigenstates are assumed to be identical with the gauge eigenstates. 



Names 


Electroweak Eigenstates 


Mass Eigenstates 


Higgs 




h, H, A°,G°,H^,G^ 




UL,UR,dL,dR 


Ui,U2,dl,d2 


squarks 


SL,Sb,,Cl,Cr 


Sl, S2, Cl, C2 




tL,tR, bL, bji 


tl, <2, &2 




CL, CR,, v e 


ei,e2,Ve 


sleptons 








TL,TR,Vr 




neutralinos 




Xl,X2,X3,X4 


charginos 






gluino 


9 


saine 



Table E.3: Gauge and Mass Eigenstates of the MSSM 
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E. 1.5 Implementation in SARAH 
E.l. 5.1 Model File 

Here, we want to explain the model file of the MSSM in SARAH. 

a) Tlie gauge sector is U{1) x SU {2) x SU{3) and is just defined by adding the corresponding vector 
superfields. 

Gauge [[1]]={B, U[l] , hypercharge, gl.False}; 
Gauge[[2]]={WB, SU[2] , left, g2,True>; 
Gauge [ [3] ] =-[G , SU[3] , color, gS.False}; 

b) The doublot suporfiolds arc q, I, Hd, Hu and addcd by 

Fields[[l]] = {{uL, dL}, 3, q, 1/6, 2, 3}; 

Fields[[2]] = {{vL, eL}, 3, 1, -1/2, 2, 1}; 

Fields[[3]] = {{HdO, Hdm} , 1, Hd, -1/2, 2, 1}; 

Fields[[4]] = {{Hup, HuO}, 1, Hu, 1/2, 2, 1}; 

c) The singlets superfields are d, u, e and added by 

Fields[[5]] = {conj[dR], 3, d, 1/3, 1, -3}; 

Fields[[6]] = {conj [uR] , 3, u, -2/3, 1, -3}; 
Fields[[7]] = {conj [eR] , 3, e, 1, 1, 1>; 

d) The superpotential is defined by 

SuperPotential = { {{1, Yu},{u,q,Hu», {{-1 , Yd} , {d, q,Hd» , 

{{-l,Ye},{e,l,Hd}>, {{l,\[Mu]>,{Hu,Hd» }; 

e) The gauge fixing terms for the unbroken gauge groups are 

DEFINITION [GaugeES] [GaugeFixing] = 
{ {Der[VWB], -1/(2 RXi[VWB])}, 

{Der[VG], -1/(2 RXi [VG] ) »; 

f) The name of the new eigenstates is 'EWSB': 
NameOfStates=-[GaugeES, EWSB}; 

g) The vector bosons and gauginos rotate as foUows 
DEFINITION [EWSB] [GaugeSector]= 

{ {VWB, {l,{VWm,l/Sqrt[2]},{conj [VWm] ,1/Sqrt[2]», 

{2 , {VWm, -\ [Imaginaryl] /Sqrt [2] } , {conj [VWm] , \ [Imaginaryl] /Sqrt [2] }} , 

{3,{VP, Sin[ThetaW]},{VZ, Cos [ThetaW] }}}, 
{VB, {1,{VP, Cos [ThetaW]},{VZ,-Sin [ThetaW]}}}, 
{fWB, {l,{fWm,l/Sqrt[2]}, {fWp, 1/Sqrt [2] }}, 

{2 , {f Wm, -\ [Imaginaryl] /Sqrt [2] } , {f Wp , \ [Imaginaryl] /Sqrt [2] }} , 

{3,{fW0,l}}} 

}; 
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h) The neutral components of the scalar Higgs get a VEV 



DEFINITION [EWSB] [VEVs] = 

{■[SHdO, {vd, 1/Sqrt[2]}, {sigmad, \ [Imaginaryl] /Sqrt [2] } ,{phid, 1/Sqrt[2]», 
{SHuO, {vu, 1/Sqrt[2]>, {sigmau, \ [Imaginaryl] /Sqrt [2] },-[phiu, 1/Sqrt[2]>»; 



i) The particles mix after EWSB to new mass eigenstates 

DEFINITION [EWSB] [MatterSector] = 
{{{SdL, SdR}, {Sd, ZD}}, 

{{SuL, SuR>, {Su, ZU}}, 

{{SeL, SeR}, {Se, ZE», 

{■Cphid, phiu}, {hh, ZH}}, 

{{sigmad, sigmau}, {Ah, ZA}}, 

{{SHdm, conj [SHup] } , {Hpm, ZP}} , 

{{fB, f WO, FHdO, FHuO}, {LO, ZN}}, 

{{{fWm, FHdm>, {fWp, FHup», {{Lm,Um}, {Lp,Up»»; 



j) The new gauge fbcing terms are 



DEFINITION [EWSB] [GaugeFixing] = 

{{Der[VP], - 1/(2 RXi[VP])}, 

{Der [VWm] +\ [Imaginaryl] Mass[VWm] RXi [VWm] Hpm[{l}] , - l/(RXi [VWm] )}, 

{Der[VZ] + Mass[VZ] RXi [VZ] Ah[{l}] , - 1/(2 RXi[VZ])>, 

{Der [VG] , - 1/ (2 RXi [VG] ) » ; 



k) No particles should be integrated out or deleted 



IntegrateOut={} ; 
DeleteParticles={} ; 



1) The Dirac spinors are the following 



dirac[[l]] = 


{Fd, 


dirac [[2]] = 


{Fe, 


dirac [[3]] = 


{Fu, 


dirac [[4]] = 


{Fv, 


dirac [[5]] = 


{Chi, 


dirac [[6]] = 


{Cha, 


dirac [[7]] = 


{Glu, 


dirac [[8]] = 


{Bino 


dirac [[9]] = 


{Wino 


dirac [[10]] = 


= {HDO 


dirac [[11]] = 


= {HUO 


dirac [[12]] = 


= {Hp, 


dirac [[13]] = 


= {Hm, 



FdL, FdR}; 

FeL, FeR}; 

FuL, FuR}; 

FvL, 0}; 

LO, conj[LO]>; 

Lm, conj [Lp] } ; 

f G, conj [fG]>; 
, f B, conj [fB]}; 
, fWB, conj [fWB]}; 
, FHdO, conj [FHdO]}; 
, FHuO, conj [FHuO]}; 

FHup , conj [FHup] } ; 

FHdm , conj [FHdm] } ; 



m) A M^file and a FeynArts model file should be written automatically for the eigenstates after 
EWSB 
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makeOutput = { {EWSB, -[TeZ, FeynArts}} }; 

n) The spectrum from the file SPSI . spc should be read in 
SpectruinFile= "SPSI. spc"; 

E. 1.5. 2 Names of Particles and Parameters 

In the following we list all names of fields appearing in the MSSM in SARAH. 
• Gauge Eigenstates 

— Fermions 



d? = 





FeL[{gen>] \ 
FeR[{gen}] J 
FvL[-[gen}] \ 

FdL[-Cgen,col}] 
FdR[{gen,col}] 
FuL [{gen, col}] 
FuR [{gen, col}] 




HDO 



HUO 



HDM = 



HUP = 



FHdO 
conj [FHdO] 

FHuO 
conj [FHuO] 

FHdm 
conj [FHdm] 

FHup 
conj [FHup] 



B = 
W' = 
G' = 



Bino 



Wino = 



Glu = 



y conj [f B] 
f W [{gen}] 
conj [f W [{gen}]] 

f G [{gen}] 
conj [f G [{gen}] ] 



— Scalars 
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SeR[{gen}] 




SeL [{gen}] 








SvL [{gen}] 




SdR[{gen,col}] 


di 


SdL[{gen,col}] 




SuR[{gen,col}] 


dT 


SuL [{gen, col}] 




SHdO 


Ha 


SHdm 


H" 


SHuO 


H+ 


SHup 



— Vector Bosons and Ghosts 



B, 


VB [{lor}] 








VW[{gen,lor}] 




gW[{gen}] 


4 


VG [{gen, lor}] 


4 


gG [{gen}] 



• Eigenstates after Electroweak Symmetry Breaking 
— Fermions 



— Scalars 



d? 




G' = 




FeL[{gen}] \ 
FeR[{gen}] J 
FvL[{gen}] \ 

O J 
FdL [{gen, col}] 

FdR[{gen,col}] 

FuL [{gen, col}] 

FuR[{gen,col}] 



Chi [{gen}] 
Cha[{gen}] 



LO [{gen}] 
conj [LO [{gen}] ] 
Lm [{gen}] 

conj [Lp [{gen}] ] 



Glu: 



f G [{gen}] 
conj [f G [{gen}]] 



e' Se [{gen}] 
d''" Sd [{gen, col}] 


i>i SvL [{gen}] 
u*'" Su [{gen, col}] 


hi hh[{gen}] 
H f Hpm[{gen}] 


Ah [{gen}] 



— Vector Bosons and Ghosts 
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7m 


VP [{lor}] 








VZ[{lor}] 


r]z 


gz 




VWm [{gen}] 


, 'nw+ 


gWm,gWinC 


s; 


VG[{gen,lor}] 




gG [{gen}] 



— Mixing Matrices 





ZD[i,j] 


yU 


ZU[i,j] 


yE 


ZE[i,j] 


yA 


ZA[i,j] 




ZH[i,j] 




ZP[i,j] 


yN 


ZN[i,j] 




Um[i,j] 




Up[i,j] 



— Couplings and other parameters 



i j 


Yd[i,j] 


Ad 


A[Yd] [i, j] 




\[Mu] 


i j 


Yu[i,j] 


AU 


A[Yu] [i, j] 


B, 


A[\[Mu]] 




Ye[i,j] 




A[Ye] [i, j] 


Qw 


ThetaW 




gl 


92 


g2 


93 


g3 


Vd 


vd 


Vu 


vu 


v 


v 



E.2 MSSM in CKM Basis 

Name in SARAH: MSSM-CKM 

The first implementation of the MSSM in SARAH deals with flavor violation only in the SUSY sector if of 
diagonal clcmonts of the soft breaking terms are not forbidden. But it is known that also in the Standard 
model sector flavor violation appears in the Yukawa sector. The convenient choice for the parametrization 
is the so called CKM basis. In this section we want to demonstrate that it is easy to incorporate flavor 
violation in the Standard model sector, i.e. we show how to deflne the MSSM in the CKM basis in SARAH. 



E.2.1 Motivation 

The Yukawa matrices are not necessarily diagonal in the gauge eigenstates. This will lead after EWSB 
to a rotation of the gauge eigenstates to a new basis, which diagonalizes the Yukawa matrices: 

dL = VJl, UL=Vuf?j, (E.32) 

4 = Udv'l, UR = Uuul, (E.33) 

where we have assigned the gauge eigenstates with a upper index 0. We can define diagonal Yukawa 
Yu, Yd matrices in this new basis 

Frf = ulY^Vd, Ya = UlY^Vd (E.34) 
In SUSY theories the soft breaking parameters are rotate in a similar way: 

= Vd^fVd, ml = UimfUu, ml = UlmfUd (E.35) 

Ad = U^A^Vd, A^ = UlAl'^Vu (E.36) 

Aftcr this redefinition the mixing matrices Vd,Vu,Ud,Uu disappear from the Lagrangian and only the 
product of Vd and Vu stays: 

VcKM = V^Vd (E.37) 
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E. 2. 2 Implementation in SARAH 
E.2.2.1 Model File 

We point out only the differences to the first model file of the MSSM: 

a) We name the gauge eigenstates with a index 0: 

Fields[[l]] = {{uLO, dLO}, 3, qO, 1/6, 2, 3}; 

Fields[[5]] = {conj [dRO] , 3, dO, 1/3, 1, -3>; 
Fields[[6]] = {conj [uRO] , 3, uO, -2/3, 1, -3}; 

b) The name of the eigenstates are 'SCKM' and 'EWSB' 
MameOf States={GaugeES , SCKM , EWSB} ; 

c) We first rotate quarks and squarks to the CKM basis. If a CalcHep model should be produced, it 
is helpful also to rotate the auxiliary fields. 

DEFINITION [SCKM] [MatterSector] = 
{{{SdLO}, {SdL, Vd}}, 
{{SuLO}, {SuL, Vu», 
{{SdRO}, {SdR, Ud}}, 
{{SuRO}, {SuR, Uu», 
{{AdLO}, {AdL, Vd», 
{{AuLO}, {AuL, Vu», 
{{AdRO}, {AdR, Ud», 
{{AuRO}, {AuR, Uu», 

{{{FdLO}, -CFdRO», {{FdL,Vd}, -[FdR,Ud»}, 
{{{FuLO}, {FuRO», {{FuL,Vu}, {FuR,Uu»»; 

d) That's all, now we can go an like in the first case 
E. 2. 2. 2 Parameter File 

In the parameter file we can do the following definition: 
a) Rotation of the Yukawa matrices 

{YuO, {Dependence -> suin[i001 , 1 ,3] *sum[i002, 1 ,3] *Yu[i001 , 1002] *Delta[i001 , 1002]* \ 

Vu [1001 , lndexl] *conj [Uu [1002 , lndex2] ] » , 



b) Rotation of the soft breaking masses 

{mq02, {Dependence -> sum[1001 , 1 ,3] *suin[1002, 1 ,3] *mq2 [1001 ,1002] *Delta[1001 , 1002]* \ 

Vd[1001,lndexl]*conj [Vd[1002,lndex2]] », 



c) Definition of the CKM matrix as product of two matrices 

{CKM, {MatrlxProduct -> {Vd,Vu}, 
...» 
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E. 3 Next to Minimal Supersymmetric Standard Model 

Name in SARAH: NMSSM, One_Rotation and NMSSM, Two_Rotations 



E. 3.1 Additional Gauge Singlet and Superpotential 

The second simplest, realistic SUSY model after the MSSM is the so called NMSSM (see ^ for a 
review). In this model a gauge singlet S is added to the particle content. This singlet couples to the 
Higgs superfields in the superpotential and also has a cubic self coupling: 



Wnmssm 



-LYMd - QYddHd + QYuuHu + \SHuHd + t:'^^'- 



(E.38) 



The last cubic term is added to forbid a Peccei-Quinn symmetry. The singlet field also receives a VEV 
because of its mixing with the Higgs fields 



v2 



(E.39) 



Obviously, the VEV and the coupling A of the superpotential plays the role of the /x-term in the MSSM 



{S) 
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(E.40) 



E.3.2 Higgs Sector 

The scalar and the pseudo scalar component of the singlet fields mixes with the Higgs field. The new 
basis in the gauge eigenstates are — (0^, i^^, (ps)^ and a — (ct^, ct^, cTs)^- But is also common to isolate 
first the Goldstone by 



(E.41) 



and rotate afterwards (cr, cts). 
E. 3. 3 Neutralino Sector 

The fermionic component of the singlet mixes after electroweak symmetry breaking with the neutral 
Higgsinos. This extends the neutralino mass matrix to five dimensions and can be written in the basis 



/ 






^ cos/S 


sin /3 






( \ 








— sin /3 


cos/3 







a 


v 


0-S j 




[ 





1 ) 




\ ) 



1 


Ml 










M2 




— ei>5 


e,v\ 




2cw 


2sw 






— eV2 




2cw 


2sw 


\ 









— evi 


eV2 


2cw 


2cw 


evq 


-2V2 




2sw 





-/^ 


-M 
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O 

o 



V2 



Vd 



(E.42) 



\/2kVs j 
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E. 3. 4 Implementation in SARAH 
E.3.4.1 Model File: One Rotation 

For the NMSSM the MSSM model file has to be changed as foUows: 

a) The gauge singlet is added to the particle content 

Fields[[8]] = {sR, 1, s, O, 1, 1}; 

b) Two new trilinear term are added to the superpotential 

{..., {{lA[Lambda]},{Hu,Hd,s}}, 

{{l/3,\[Kappa]},{s,s,s}} }; 

c) The gauge singlet receives also a VEV 

{SsR, {vS, 1/Sqrt[2]}, {sigmaS, \ [Imaginaryl] /Sqrt [2] }, 
{phiS, 1/Sqrt[2]» } 

d) There are new mixing between the gauge singlet and the other particles 

{{phiu, phid, phiS}, {hh, ZH}}, 

{{sigmau, sigmad}, {AhT, ZT}}, 

{{fB, f WO, FHdO, FHuO.FsR}, {LO, ZN}}, 

E.3.4.2 Model File: Two Rotations 

It is also common to define two rotations in the pseudo scalar sector: First rotate the pseudo scalar Higgs 
to the Goldstone before doing the mixing with the pseudo scalar of the singlet. This can be done by 

a) Two rotations are initialised 

NameOf States=-[GaugeES, TEMP,EWSB> 

b) First a temporary rotation 

DEFINITION [TEMP] [MatterSector] = 
{{sigmau, sigmad}, {AhT, ZT»; 

c) Afterwards the rotation to the EWSB eigenstates 

DEFINITION [EWSB] [MatterSector] = 
{{AhT, SigmaS}, {Ah, ZA}} 

d) Proper definition of the mixing matrices: 

{ZT,{Dependence -> {{Cos [\ [Beta] ] , Sin[\[beta]]}, 

{-Sin [\ [Beta] ] , Cos [\ [beta] ] }}} , 

{ZA, {Dependence -> {{1, O, O }, 

{O, Cos[\[Phi]], -Sin[\[Phi]]} 
{O, sin[\[Phi]], Cos[\[Phi]]}}} 
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E. 3. 4. 3 Names of Particles and Parameters 

There are only small changes in the particle definitions of the NMSSM in comparison to the MSSM file. 
The names of the mass and are the same, just to the gauge eigenstates the singlet is added: 



S Ss 



Sing = 



E. 4 MSSM with bilinear R-parity Violation 

Name in SARAH: MSSM-BiRpV 
E.4.1 Introduction 

R-parity was introduced in SUSY to forbid several couplings and protect the proton from decay. But, it 
is possible to allow some coupUngs which brcak R-parity without running into trouble. For example, the 
Higgs Hd and the lepton doublet have the same gauge quantmn numbers, and therefore the bihnear term 

Wrpv = eiHji (E.43) 
can be added to the superpotential. 
To the Lagrangian also the soft breaking couphngs 

B,,Hji (E.44) 
and a mass term 

mlH^,H*Ji + c.c. (E.45) 
are added. 

Because of their mixing with the Higgs the sneutrinos will receive a VEV after EWSB 

= (0i + icr^ + v^) (E.46) 

E.4.2 Neutrino and Neutralino Sector 

A nice aspcct of bilinear R-parity violation is. that it Icads to a mixing between the neutrinos and neutrali- 
nos and it offers a possible explanation for neutrino masses. The new basis is = {B, H^, v^, t'/j^'r) 
and leads to seven max eigenstates Ao: 

Xi = Z,''ji.^ (E.47) 
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Z N diagonalizes the combined neutriiio/neutraliiio mass inatrix 
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E. 4. 3 Chargino and Electrons 

Because of the bilinear R-parity violation, also the charginos do mix with the charged leptons. The new 
basis is = {HJ, W~ , e^, fi^, tl) and , cr, fiR, tr) and we have the followmg relations: 



(E.49) 
(E.50) 
(E.51) 



With the mixing matrices Z and Z^ can the mass matrix M be diagonahzed 
/ M 



M = 



9'"u 
72 

O 

o 
o 



92 f d 
V2 

M 

Y^^V2 

yr^>3 

V2 



92 f 1 
V2 

-£l 

-O 

-o 



92 "2 

72 

-£2 
O 

V2 

O 



92^3 \ 

n/2 * 

-£3 

O 

o 

V2 / 



(E.52) 



E. 4. 4 Higgs and Sneutrino/Selectrons 

In the scalar sector we have the new mixings 

• Charged scalars 

The basis for charged scalars is — {H^ * , 6^,6/?) and leads to a 8 x 8 mass matrix. 

• Neutral scalars 

The basis for neutral scalars is — Re {H^j, H^, Dl) 

• Pseudo scalars 

The basis for neutral pseudo scalars is P" = Jm [H^, i'l) 

The corresponding mass matrices can be found e. g. in |27) . 



E. 4. 5 Implementation in S ARAH 
E.4.5.1 Model File 

To add bilinear R-parity violation to the MSSM, the following changes are needed 
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E.4. MSSM WITH BILINEAR R-PARITY VIOLATION 



a) Add the corresponding term to the superpotential 
SuperPotential={ . ■ ■ , {{l,\[Epsilon]},-[l,Hu}}}; 

b) Add the non canonic soft breaking mass 

DEFINITIONCGaugeES] [Additional] = {{{conj [SvL] , SHdO}, {1, mHL2}}, 
{{conjESeL], SHdm}, {1, mHL2»>; 

c) Give the sneutrinos a VEV 

{SvL, {vL, 1/Sqrt[2]}, {sigmaL, \[ImaginaryI]/Sqrt[2]}, 
{phiL, 1/Sqrt[2]}} >, 

d) Change the mixings 



{{phiu, phid,phiL}, {hh, ZH}}, 
{{sigmau, s igmad, sigmaL}, {Ah, ZA}}, 
{{conj [SHup] ,SHdin,SeL,SeR},{Hpin,ZP}}, 



E.4.5.2 Names of Particles 

The mass eigenstates in the MSSM model with R-parity violation are 
• Fermions 



• Scalars 



d? = 



G' = 



'l, a 
L 

■i, a 
R 
i,a 
L 
i, a 



A- 



'h 

2il 



Fd[{gen,col}] = 
Fu[{gen,col}] = 



FdL [{gen, col}] 
FdR[{gen,col}] 
FuL[{gen,col}] 
FuR[{gen,col}] 



Chi [{gen}] = 
Cha[{gen}] = 



LO[{gen}] 
conj [LO [{gen}] ] 

Lm [{gen}] 
conj [Lp [{gen}] ] 



Glu = 



f G [{gen}] 
conj [f G [{gen}]] 



Sd [{gen, col}] 


u''°' Su [{gen, col}] 


hi hh [{gen}] 
Hf Hpm[{gen}] 


A'y Ah [{gen}] 



Vector Bosons and Ghosts 
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lu 


VP [{lor}] 


rjp 




z. 


VZCUor}] 


1)2 


gz 




VWm[{gen}] 


VW- , Vw+ 


gWm,gWmC 


5; 


VG [{gen, lor}] 


4 


gG [{gen}] 



E.5 //i/SSM 

Name in SARAH: munuSSM 

E.5.1 Superpotential 

The so called fjbvSSM combines the idea of the NMSSM with bilinear R-parity violation: The /x-term 

is generated dynamically by the VEV of a gauge singlet and neutrino masses can be explained because 
of the niixing with neutrahnos. The gauge singlet in the /zi^SSM is the right handed neutrino. The 
superpotential is 



- - I ^ 

W = -LYeeHd - QYddHd + QYy,uHu + LY^Hu^r - XvRHuHd + -«£-1 



(E.53) 



The number of generations of vr varies, 1 or 3 are used. We consider in the following only the cases with 
one generation. 

The scalar gauge singlet Dr and the left handed sneutrinos mix with the Higgs and receive also a VEV 
after EWSB: 



i^R = —?E i(l>R + i<^R + Vr) 



V2 
1 

71 



{(pi + icji + Vi) 



This leads to the following definition of the /U term 

M = ^XVR 



(E.54) 
(E.55) 

(E.56) 



E. 5. 2 Neutralino and Neutrino Sector 

The superpotential of the /ui^MSSM breaks R-parity and leads to a mixing between the neutralinos and the 
neutrinos. Therefore the basis of neutral fermions \s (f) = {B, W^, H^, H^, ur, vl) and the corresponding 
mass matrix has the form 



M: 



m 



X,VR 

T 



M, 



m 



D 



O 



(E.57) 



Here, is the mass matrix of the the neutralinos in the MSSM, m-^^^VR is the mixing between neutralino 
and gauge singlet 



0,0, -^\vu,-^\vd, -^^WR 



(E.58) 
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E.5. fiiySSM 



is the mixing between neutralinos and neutrinos 



O £1 ^ 



(E.59) 



with 



The Dirac term is 



(E.60) 



m'o = ^i">« (E.61) 
E. 5. 3 Chargino and Electron Sector 

The charged leptons niixes with charged compoiients of the Higgsinos and gauginos, so the basis is 
if)- = [W~ ,H'j^ jCl) and ijj^ = {W^ ,ej^). This leads to a 5 x 5 mass matrix 
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(E.62) 



E. 5. 4 Higgs and Selectron/Sneutrinos 

Also the Higgs and the selectron/sneutrino sector mixes. This leads to the fohowing rotations: 

• Charged scalars 

The basis for the mixing of the charged scalars is — {HJ'* , ,eL,eji) and leads to a 8 x 8 
mass matrix. 

• Neutral scalars 

The basis for the mixing of the neutral scalars is S'" — Re (-ff", i?,", i'r, vl) 

• Pseudo scalars 

The basis for mixing of the neutral pseudo scalars is P'^ — Im {H^, iJ^, i/R, ul) 
AU corresponding mass niatrices can be found in |28| . 

E. 5. 5 Implementation in SARAH 
E.5.5.1 Model File 

The model file for the /ii/SSM can be built by using the model file for the MSSM and making the following 
changes 



a) Add the right handed neutrino superfield 
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Fields[[8]] = {vR, 1, v, O, 1, 1}; 

If you use conj [vR] or vR depends on your coiivention of the assigned lepton number to vR. 

b) Add three terms to the trilinear superpotential 

SuperPotential = {...{1, Hu, v, Yv, 1}, {v ,v, v, \ [Kappa] , 1/3} , {Hu, Hd, v, \ [Lsunbda] , -1}} ; 

c) Give the right and left haiided sneutrinos a VEV 

{SvL, {vL, 1/Sqrt[2]}, {sigmaL, \ [Imaginaryl] /Sqrt [2] } , 

{phiL.l/ Sqrt[2]}} } 
{SvR, {vR, 1/Sqrt[2]}, {sigmaR, \ [Imaginaryl] /Sqrt [2] } , 

{phiR.l/ Sqrt[2]» } 

d) Change the mixings 

{{phid, phiu, phiR, phiL} , {hh, ZH}}, 
{{sigmad, sigmau, sigmaR, sigmaL}, {Ah, ZA}}, 
{{SHdm, conj [SHup] , SeL, SeR} , {Hpm, ZP}} , 
{{fB, f WO, FHdO, FHuO, FvR, FvL} , {LO, ZN}}, 

{{{fWm, FHdm, FeL} , {fWp, FHup, conj[FeR]}}, {{Lm,Um},{Lp,Up}}} 
E. 5. 5. 2 Names of Particles 

The right handed neutrino is added and the definition of the Dirac spinor for neutrinos is changed. 

FvL [{gen}] \ 
conj [FvL [{gen}] ] j 

The names of the eigenstates after electroweak symmetry breaking are the same as in the oase of the 
MSSM with bihnear R-Parity violation. 



vr SvR 



Fv[{gen}] 



E.6 Other Models 



We have also added some more models to the package in order to give examples for different possibihties. 



E.6.1 MSSM with Seesaw II 
E.6. 1.1 Motivation 

For checking the RGEs of models including non-fundamental representations, we have used a SU{^)- 
inspired extension of the MSSM [221 [30]- An additional 15 and 15 of SU{5) is added. The 15 decomposes 
in 

15^S + T + Z (E.63) 



84 



E.6. OTHER MODELS 



The fields have the charges 

S: (m, 4; 

T: (1,3,1) 
1 



Z : ^3,2, g 
uiider SU{3) x SU{2) x U{1). 

The superpotential terms after 5t/(5)-breaking are 

Wmssm + (YtLTL + YjSd) + yJzL 



V2 

+ -i= (\iHdfHd + \2HufHu) + MtTT + MsSB + MzZZ 

Integrating out the new, heavy states leads to efFective operators for the neutrino masses: 
1 



Yt\2vi (E.64) 



2Mt 

Remark Our one loop RGEs are in agreement to [30], but not to |29| . 

E. 6. 1.2 Implementation in SARAH 

a) The new fields are 

Fields[[8]] = {{{T0,-1/Sqrt[2] Tp} , {-l/Sgrt [2] Tp, -Tpp}}, 1, t, 1, 3, 1}; 
Fields[[9]] = {{{Tmm, -1/Sqrt [2] Tm} , {-1/Sqrt [2] Tm, -TOb», 1, tb, -1, 3, 1>; 
Fields [[10]] = {S, 1, s, -2/3, 1,6}; 
Fields[[ll]] = {conj [Sc] , 1, sb, 2/3, 1, -6}; 
Fields[[12]] = ■[■[zl,z2}, 1, z, 1/6, 2, 3}; 
Fields[[13]] = {{zlb,z2b}, 1, zb, -1/6, 2, -3}; 

b) The superpotential is 

{ {{1, Yu},{u,q,Hu}}, {{-l,Yd},{d,q,Hd», 
{{-l,Ye},{e,l,Hd», {{1,\ [Mu] } , {Hu , Hd}} , 
{{1/Sqrt[2] ,Yt},{l,t,l}}, {{1/Sqrt[2] , Ys} , {d, s , d}} , 
■[■[l,Yz},-[d,z,l}}, {■[1/Sqrt[2] ,Ll},{Hd,t,Hd}}, 
{■[1/Sqrt[2] ,L2},{Hu,tb,Hu», {{I ,MT} , {t , tb» , 
{{1 , MZ} , -Cz , zb}} , {{1 , MS} , {s , sb}}} ; 

c) The new field are integrated out 

IntegrateOut= {Szlb , Sz2b , Szl , Sz2 , SSc , SS , STO , STmm , 
Fzlb , Fz2b , Fzl , Fz2 , FSc , FS , FTO , FTmm} ; 

E.6.2 Standard Model 

Name in SARAH: SM 
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E.6.2.1 Motivation 

As example, how to get a non supersymmetric limit of a SUSY theory using SARAH, we demonstrate, 
how to get the Standard model. This procedure can surely be used for other models, like double Higgs 
model. 

E.6.2.2 Implementations in SARAH 

Model File To get the SM, do the foUowing with the MSSM 

a) First, all non Standard model Particles have to be deleted: 

DeleteParticles = {SciR, SdL, SuR, SuL, SeR, SeL, SvL, fB, fW, fG, hh[{2,2}] , Ah, Hpm}; 

b) As second step, the Higgs self coupling must be changed 

DEFINITION [EWSB] [Additional] = {{hh,hh,hh,hh} , {1/24, \lambda}}; 

Names of Fields in SARAH The names are the same as for the MSSM, just without all SUSY 
particles. 

E.6.3 Effective Theory: MSSM with very heavy Gluino 
E.6.3.1 Motivation 

In particular SUSY scenarios the gluino is very heavy and it can be a appropriate choice to consider only 
a effective theory by integrating out the gluino. 

E. 6. 3. 2 Implementation in SARAH 

Model File For integrating out the gluino, you have just to add 

IntegrateOut[[l]] = f G 

to the model file of the MSSM. 

Names of Particles The names of particles are exactly the same as for the full MSSM. 

E. 6. 4 Additional [/(l)-Gauge Group 
Name in SARAH: MSSM-Ul 

E.6.4.1 Motivation 

This is a extension of the MSSM which should demonstrate, how easy it is to add a new gauge interaction 

in SARAH. Thorcforc, wc want to extend the MSSM by a new abelian gaugo group, just called extra. 
The coupling constant should be g4 and only the leptons are charged under this gauge group. 
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E.6.4.2 Implementation in SARAH 



The Model File To build a MSSM with an additional gauge group do the foUowing 



a) Add the vector superfield 



Gauge [ [4] ] ={X , 1, extra, 



g4, False}; 



b) Give your particles a new quantum number 

Fields[[l]] = {{uL, dL} , 3, q, 1/6, 2, 3, 0}; 

Fields[[2]] = {{vL, eL}, 3, 1, -1/2, 2, 1, l/2>; 

Fields[[3]] = -C-CKdO, Hdm}, 1, Hd, -1/2, 2, 1, 1/2}; 

Fields[[4]] = {-CKup, HuO}, 1, Hu, 1/2, 2, 1, -1/2}; 

Fields[[5]] = {conj [dR] , 3, d, 1/3, 1, -3, 0}; 
Fields[[6]] = {conj [uR] , 3, u, -2/3, 1, -3, 0}; 
Fields[[7]] = {conj [eR] , 3, e, 1, 1, 1, -1/2}; 

This is only a toy model to clarify the possibility of defining new gauge groups in SARAH. This 
model is not free of gauge anomalies! This can also be seen by all the warnings SARAH prints by 
evaluating this model. 

c) The new gaugino mixes with the neutralinos because the Higgs are charged under the new gauge 
group 

{{fB, f WO, FHdO, FHuO, fX}, {LO, ZN}}, 

d) Define a new majorana fermion for the gauge eigenstates 
dirac[[14]] = {Xino, fX, conj [f X]}; 

Names of Particles The names are the same as for the MSSM plus a new vector boson and gaugino 
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Appendix F 

Verification of Output 



We have tested the model files generated by SARAH for the MSSM after electroweak symmetry with the 
existing model files of FeynArts and CalcHep/CompHep in several ways. 

F.l FeynArts 

F. 1.1 Comparison on Vertex Level 

It is quite easy to translate the Standard MSSM model file included in FeynArts to the parametrization 
used by SARAH by because both are packages for Mathematica. So our first check was comparing the 
vertices by doing the foUowing steps: 

a) Generate a model file from SARAH 

b) Read the model file and the Standard FeynArts model file 

c) Expand all indices (leads to more than 15000 interaction!) 

d) Replace the parameters from FeynArts with the corresponding parameters of SARAH 

e) Insert numcrical values (of SPS 1) 

f) Check that the same interaction have the same numerical values at least up to eight digits. 

F. 1.2 Compciring on Amplitude SquaLred Level 

The second step of checking our FeynArts model file was to calculate the amplitude squared of several 
proccsscs to makc sure that there is no problem with intcrfcrcnccs. Hcrc wc have chcckcd about 50 
processes (more have been done with CalcHep) with different generic structure and compared again the 
numerical values. For 2^2 processes we have set s = (lOOOOGeF)^ and t = (lllGeV)^ to get a 
numerical values for the amplitude squared. 

F. 2 CalcHep/MicrOmegas 

With CalcHep and MicrOmcgas wo did tho thrco stcps to verify our produced model file. But first, we 
must say some words about the used model file of CalcHep. 
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F.3. RENORMALIZATION GROUP EQUATIONS 



F.2.1 Used Model File for Checks 

In the MSSM file used by the current version of CalcHep the one loop corrections for the Higgs self 
interactions are added. Of course, this can't be done by SARAH. Furthermore, the interactions of 
the Goldstone bosons are related to the interactioii of the corresponding gauge boson and aren't the 
interaction of the pseudo scalar. 

In order to be comparable with the model file produced by SARAH we have used a older version of a 
MSSM file of CompHep, which can be found under theory. sinp.msu.ru/~semenov/inssm.html. In this 
file it is possible to switch the loop corrections ofF and the Goldstone interaction are the interaction of 
the first pseudo scalar Higgs. 

F. 2. 2 Comparison on Vertex Level 

For comparing the model files on vertex level, we first formatted them in a way to be readable in 
Mathematica by using a SED script. After that we used relations between the parameters of CalcHep 
and the parameters of SARAH and introduces numerical values for the free parameters. One must 
be aware of the fact that the convention for the mixing matrices in SARAH and CalcHep differs by a 
transposition. 

F. 2. 3 Comparison on Process Level 

The second step in the validation of our CalcHep model was to calculate more than 200 cross sections and 
decay widths using the Vegas routine of CalcHep. About 40 calculations have been made in unitarity- 
as well as in Feynman-gauge. In addition, more than 20 four point interactions generated by auxiliary 
fields in CalcHep and SARAH have been separately checked by switching off the other diagrams. 

F. 2. 4 Comparison of Relic Density 

The last step was calculating the relic density with MicrOmegas for SPS 1 and SPS 2. Besides the values 
for flh^ the contribution of the different annihilation channels was checked: The maximal difference for 
all channels contributing more than 0.01 % was between the dominating channel of SPS 2 (xiXi ~^ h^)'- 
30.00% to 29.99%. 

F. 3 Renormalization Group Equations 

We have compared the analytical results for the one and two loop RGEs calculated by SARAH for the 
MSSM with |[3T], for the NMSSM with p2] and for the MSSM with bihnear R-Parity violation with [33] . 
In addition, we have checked a model containing non-fundamental representations: the SU{5) inspired 
Seesaw II model of [30] and [53]. It is know that the there are discrepencies of the RGEs given in this 
two papers. Our result fully agrees with |30j . 

Also numerically checks have been done by comparing the RGEs for the MSSM with the RGEs imple- 
mented in SPheno [M]. Both set of RGEs were in fuU agreement. 



89 



Appendix G 

Evalution Time 



To give a impression for thc needed evaluation time for different routines and models, we collected some 
values iii table IG.ll This times were measured under Mathematica 5.2 running oii a Intel Quad CPU 
Q8200 with 2.33 GHz and 4GB RAM. 



Command 


MSSM (No FV) 


MSSM-CKM 


NMSSM 1 


fivSSM 


MSSM + [/(1) 


Start 


12.75 


18.03 


19.02 


27.06 


16.14 


ModelOutput [EWSB] 


74.83 


78.70 


94.64 


115.08 


110.47 


MakeFeynArts 


0.74 


3.58 


1.12 


0.98 


0.48 


MakeCalcHep [] 


6.03 


22.74 


15.57 


47.08 


6.7 


MakeTeX [] 


0.81 


5.79 


1.25 


1.38 


1.48 


CalcRGEs [] 


50.72 


50.8 


91.07 


265.29 


68.18 


CalcLoopCorrections [EWSB] 


7.07 


28.44 


8.14 


7.98 


8.84 


Vertex[{bar [Fd] ,Fd,VP}] 


0.05 


0.06 


0.06 


0.06 


0.05 


Vertex [{Sd , Sd* , Su , Su*}] 


0.54 


0.61 


0.56 


0.56 


0.52 


Vertex [{VG , VG , VG , VG>] 


0.03 


0.03 


0.04 


0.03 


0.03 


Vertex[{Se,Se,VZ}] 


0.37 


0.37 


0.45 


1.13 


0.37 


Vertex [{bar [gWmC] ,gP,VWm*}] 


0.02 


0.02 


0.02 


0.02 


0.02 


Vert ex [{Hpm , VWm* , VZ}] 


0.16 


0.17 


0.26 


0.23 


0.21 


TreeMass[VZ,EWSB] 


0.15 


0.19 


0.20 


0.25 


0.48 


TreeMass [Su.EWSB] 


0.25 


0.27 


0.29 


0.34 


0.23 



Table G.l: Time needed in seconds to evaluate several commands of SARAH in Mathematica 5.2 



^With two rotations in the pseudo scalar sector 
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